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ABSTRACT 

This  report  contains  user  instructions  and  program  listings 
for  a  Loran-C  position  determination  algorithm  for  use  on  an 
Hewlett-Packard  HP-75C,  a  Sharp  PC-1500  (Radio  Shack  TRS-80  PC-2) 
and  a  Radio  Shack  Model  100  portable  computers. 


I.  INTRODUCTION 

This  report  contains  a  Loran-C  position  determination  program 
for  each  of  the  following  portable  computers:  the  Hewlett-Packard 
HP-75C,  the  Sharp  PC-1500  (Radio  Shack  TRS-80  PC-2)  and  the  Radio 
Shack)  TRS-80  Model  100.  These  programs  are  part  of  a  set  that  has 
been  developed  for  a  portable  computer  evaluation. 

The  Loran  system  is  a  radio  aid  to  navigation  which  utilizes 
the  principle  of  hyperbolic  fixing.  The  locus  of  points  for  which 
the  difference  in  arrival  time  of  synchronized  signals  from  a  pair 
of  transmitters  is  constant  determines  a  hyperbolic  line  of 
positions.  The  intersection  of  two  hyperbolic  lines  of  position 
from  two  pairs  of  stations  determines  a  hyperbolic  fix.  That  two 
pairs  of  stations  are  required  for  a  fix  does  not  necessarily  mean 
that  there  are  four  separate  stations,  for  one  station  of  one  pair 
may  be  colocated  with  one  station  of  the  other  pair  forming  a 
Loran  triplet.  Triplets  may  be  joined  "end-to-end"  by  station 
colocation  to  form  a  Loran  chain.  Loran  chains  are  common  on  both 
the  East  and  West  coasts  of  the  North  American  continent. 

The  present  day  Loran-C  operates  at  100-kHz  and  is  in  use  in 
the  Atlantic,  Pacific  and  Mediterranean  areas.  The  computational 
programs  described  herein  can  be  used  for  position  determination 
with  Loran-C  triplets.  Further  information  on  the  history, 
development  and  operation  of  the  Loran  systems  may  be  found  in 
References  1  and  2.  Details  of  the  fixing  algorithm  can  be  found 
in  Reference  3. 
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In  the  programs  presented  here,  there  are  ten  user  options. 
The  computer  display  abbreviations  for  the  option  names  and  their 
definitions  are  listed  below. 

1.  LORAN  FIX.   This  fixing  routine  is  the  main  program  for 
calculating  a  LORAN-C  fix  from  indicated  time  delays  (ITDs). 

2.  ALTERNATE  SOLN.   The  alternate  solution  routine  will  allow 
the  second  LORAN-C  fix  solution  to  be  computed.  This  routine 
toggles  an  alternate  solution  flag  so  that  on  subsequent 
fixes  Option  1  will  calculate  the  alternate  solution. 

3.  INPUT  PEST  LAT/LONG.   This  option  is  used  to  store  the 
latitude  and  longitude  of  a  destination. 

4.  HEAD  &  DIST  TO  PEST.   This  option  is  used  to  compute  the 
heading  and  distance  from  the  currect  fix  to  the  destination 
stored  by  Option  3. 

5.  PREPICT  ITPs.   This  option  may  be  used  to  predict  the 
station  ITPs  that  will  be  received  at  a  given  latitude  and 
longitude. 

6.  TWO  POINT  HEAP  &  PIST.   This  option  is  similiar  to  Option  4 
except  that  it  may  be  used  to  compute  the  heading  and 
distance  from  any  origin  to  any  destination. 

7.  CALIBRATE.   Given  the  latitude  and  longitude  of  a  benchmark 
position  and  the  indicated  time  delays  from  a  LORAN-C 
triplet  received  at  the  benchmark  position,  the  stored 
station  data  are  modified  so  that  the  fix  routine  (or 
alternate  solution  routine)  will  compute  the  known  position 
from  the  same  time  delays. 
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8.  NEW  STA.  IDs.   This  option  is  used  to  load  the  two  station 
pairs  being  received  by  requesting  the  group  repetition 
interval  (GRI)  identifier. 

9.  ASF  Corrections.   When  a  portion  of  the  propagation  path 
lies  over  land,  the  additional  secondary  factor  (ASF)  should 
be  included  in  the  observed  time  delay  if  precision 
navigation  is  required.  This  option  allows  the  user  to  input 
the  ASF  correction  [Ref.  4,  5]. 

A.  S_T_Q£.   This  option  will  exit  from  the  LORAN-C  routine. 

An  additional  EELR   option  (H)  is  available  in  the  Hewlett- 
Packard  HP-75C  and  Sharp  PC-1500  (TRS-80  PC-2)  computers.  This 
option  will  review  the  ten  options  in  the  display.  The  Radio  Shack 
TRS-80  Model  100  does  not  have  this  option  because  all  ten  options 
can  be  displayed  simultaneously. 
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II.  THE  HEWLETT-PACKARD  HP-75C. 

HP-75C  USER  INSTRUCTIONS. 


The  options  are  labeled  1  through  9,  A  or  H.  The  H  (Help)  option 
will  review  options  1  -  9  &  A.  Option  6  (TWO  POINT  HEAD  &  DIST)  is 
a  stand-alone  option  which  may  be  used  by  itself.  For  LORAN 
fixing,  Option  8  must  be  selected  first  to  load  the  station  GRl's. 
If  Option  8  is  not  used,  then  any  option  requiring  the  GRI ' s  will 
automatically  enter  Option  8  first.  Option  1  (LORAN  FIX)  must  be 
selected  at  least  once  before  Option  2  (ALTERNATE  SOLN)  can  be 
selected.  Once  the  GRI ' s  are  loaded,  any  of  the  remaining  options 
may  be  used  at  any  time.  When  first  run,  the  ASF  Correction 
Factors  are  set  to  zero;  if  they  are  subsequently  entered  as 
non-zero  values  they  will  remain  set  until  either  changed 
manually,  by  selection  of  Option  9,  or  until  the  program  is  RUN 
once  more.  In  all  cases,  South  latitudes  and  East  longitudes  must 
be  negative  for  input  and  will  be  labeled  by  a  minus  sign  '-'  on 
output.  For  ease  of  reference  the  options  are  listed  by  number 
rather  than  by  sequence  of  selection. 


INSTRUCTION 


DISPLAY 


INPUT 


Run  Program 


1.  LORAN  FIX 

ITD  for  1st  GRI 
ITD  for  2nd  GRI 
Lat  of  fix 
Long  of  fix 


2.  ALTERNATE  SOLN 
Lat  of  fix 
Long  of  fix 


3.  INPUT  DEST  LAT/LONG 


NEXT  OPTION  (Help) ? 

ITD  FOR  gri  1? 
ITD  FOR  gri  2? 
LAT  =  dd°mmlss" 
LONG  =  ddd°mm,ss" 
NEXT  OPTION  (Help) ? 


LAT  =  dd°mmIss" 
LONG  =  ddd°mmlss" 
NEXT  OPTION  (Help) ? 


DEST  LAT  dd.mmss? 
DEST  LONG  ddd.mmss? 
NEXT  OPTION  (Help) ? 

4.  HEAD  &  DIST  TO  DEST 

Heading  to  dest.     HEADING  =  dd°mm,ss" 
Dist.  to  dest.      DISTANCE  =  mmm.ff  n.mi. 

NEXT  OPTION  (Help) ? 


Run  "LORAN" 


1st  ITD  [RTN] 
2nd  ITD  [RTN] 

[RTN] 

[RTN] 


2 
[RTN] 
[RTN] 


Lat  of  dest. 

[RTN] 
Long  of  dest 

[RTN] 


4 
[RTN] 
[RTN] 
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HP-75C  USER  INSTRUCTIONS  (cont.) 


INSTRUCTION 
PREDICT  ITDs 

1st  ITD 
2nd  ITD 


DISPLAY 


LAT  dd.mmss? 
LONG  ddd.mmss? 
ITD  FOR  gri  1  = 

XXXXX.XX  MS 

ITD  FOR  gri  2  = 

xxxxx.xx  ys 
Eew  Lat/Lon  or 
Eeturn  to  menu 


NEXT  OPTION  (Help) ? 

TWO  POINT  HEAD  &  DIST 

ORIGIN  LAT  dd.mmss? 

ORIGIN  LONG  ddd.mmss? 

DEST  LAT  dd.mmss? 

DEST  LONG  ddd.mmss? 


INPUT 


N 

R 


Heading  to  dest. 
Dist.  to  dest. 


CALIBRATE 


8.  NEW  STA.  IDs 


9.  ASF  Corrections 


HEADING  =  dd0mm'ss" 
DISTANCE  =  mmm.ff  n.mi. 
NEXT  OPTION  (Help) ? 


LAT  dd.mmss? 
LONG  ddd.mmss? 
ITD  FOR  gri  1? 
ITD  FOR  gri  2? 
NEXT  OPTION  (Help) ? 


1st  GRI? 
2nd  GRI? 
NEXT  OPTION  (Help) ? 


ASF  FOR  gri  1? 
ASF  FOR  gri  2? 
NEXT  OPTION  (Help) ? 


Lat   [RTN] 
Long  [RTN] 

[RTN] 

[RTN] 

repeats  Option  5 
goes  to  next  line 


Lat  of  origin 

[RTN] 
Long  of  origin 

[RTN] 
Lat  of  dest. 

[RTN] 
Long  of  dest. 

[RTN] 

[RTN] 

[RTN] 


Lat  of  fix   [RTN] 
Long  of  fix  [RTN] 
1st  ITD  [RTN] 
2nd  ITD  [RTN] 


8 
Sta.  ID  [RTN] 
Sta.  ID  [RTN] 


1st  ASF  [RTN] 
2nd  ASF  [RTN] 


A.  STOP 
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HP-75C  SAMPLE  PROBLEM  1 


You  are  sailing  off  of  the  coast  of  California.  Your 
destination  is  Moss  Landing  at  about  36°48'N  and  121°47'W.  Your 
present  ITD  from  station  pair  GRI  9940W  is  16019  and  your  ITD  from 
station  pair  GRI  9940Y  is  42585.  What  is  your  current  position  and 
what  is  the  heading  and  distance  from  your  current  position  to 
Moss  Landing?  We  assume  that  the  GRI ' s  9940W  and  9940Y  have  not 
yet  been  loaded. 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


NEXT  OPTION  (Help)? 
DEST  LAT  dd.mmss? 
DEST  LONG  ddd.mmss? 


36.48  [RTN] 
121.47  [RTN] 


Input  destination 
36°48' 
121°47' 


NEXT  OPTION 
1st  GRI? 
2nd  GRI? 


(Help)? 


8 
9940W  [RTN] 
9940Y  [RTN] 


Input  station  ID's, 


NEXT  OPTION  (Help) ? 
ITD  FOR  9940W? 
ITD  FOR  9940Y? 
LAT  =  39°14'19" 
LONG  =  115°50I52" 


1 
16019 
42585 
[RTN] 
[RTN] 


Compute  a  fix. 


Latitude  of  fix. 
Longitude  of  fix 


This  fix  is  in  Nevada.  There  are  always  two 
solutions  to  the  LORAN  equations.  The  navigator 
must  select  the  correct  one. 


NEXT  OPTION  (Help)? 
LAT  =  35° 00' 01" 
LONG  =  125° 00  '09M 


2 

[RTN] 
[RTN] 


Select  alternate  soln, 
Latitude  of  fix. 
Longitude  of  fix. 


This  fix  is  in  the  Pacific.  The  navigator 
knows  it  is  the  correct  fix  because  of  his 
estimate  of  where  he  should  be. 


NEXT  OPTION  (Help)?  4 

HEADING  =  54° 34' 11"         [RTN] 
DISTANCE  =  190.3  8  n.mi.     [RTN] 


Find  heading  &  distance 
to  destination. 
Heading  to  Moss  Landing. 
Distance  to  Moss  Landing 


NEXT  OPTION  (Help)? 


Repeat  with  new  option 
or  press  A  to  quit. 
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HP-75C  SAMPLE  PROBLEM  2 


How  far,  and  in  what  direction,  is  Corvallis,  Oregon 
°34'N,  123°16'W)  from  Cupertino,  California  (37°19'N, 


(44 

122" 02 'W)? 


DISPLAY  CONTENTS 

NEXT  OPTION  (Help)? 

ORIGIN  LAT  dd.mmss? 
ORIGIN  LONG  ddd.mmss? 
DEST  LAT  dd.mmss? 
DEST  LONG  ddd.mmss? 
HEADING  =  353  02' 59" 
DISTANCE  =  43  8.32  n.mi 


USER  RESPONSE 

COMMENTS 

6 

Two 

point  distance 

and 

heading. 

37.19  [RTN] 

37°19' 

122.02  [RTN] 

122°02' 

44.34  [RTN] 

44°34' 

123.16  [RTN] 

123°16' 

[RTN] 

[RTN] 

NEXT  OPTION  (Help) ? 


Repeat  with  new  option 
or  press  A  to  quit. 


HP-75C  SAMPLE  PROBLEM  3 


Suppose  that  you 
receive  at  35°  North, 
from  9940W  and  9940Y. 


want  to  know  what  ITD's 
125°  West  and  at  36° 27  ' 
To  determine  these  ITD 


you  would  expect  to 

North,  126° 54'  West 

' s,  proceed  as  follows 


(assume  that  the  GRI ' s  have  been  loaded  as  in  Sample  Problem  1) 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


NEXT  OPTION  (Help) ? 

5 

LAT  dd.mmss? 

35  [RTN] 

LONG  ddd.mmss? 

125  [RTN] 

ITD  FOR  9940W  =  16019.35 

us 

[RTN] 

ITD  FOR  9940 Y  =  425  84.71 

us 

[RTN] 

Hew  Lat/Lon  or 

Eeturn  to  menu? 

N 

LAT  dd.mmss? 

36.27  [RTN] 

LONG  ddd.mmss? 

126.54  [RTN] 

ITD  FOR  9940W  =  15572.32 

ys 

[RTN] 

ITD  FOR  9940Y  =  43006.15 

ys 

[RTN] 

Hew  Lat/Lon  or 

Eeturn  to  menu? 

P 

Predict 
35° 
125° 


ITD's, 


36°  27' 
126° 54' 


NEXT  OPTION  (Help) ? 


Repeat  with  new  option 
or  press  A  to  quit. 
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HP-75C  SAMPLE  PROBLEM  4 


Suppose  you  are  receiving  an  ITD  of  16308  from  9940W  and 
42  800  from  9940Y.  These  ITD ' s  would  tell  you  that  your  location  is 
36 ° 47' 55 "N  and  121° 47  *  11 "W.  However,  you  know  that  your  position 
is  benchmarked  to  be  at  36047'36"N  and  121°46 ' 58"W,  and  you  wish 
to  calibrate  your  HP-75C  so  that  the  ITD's  of  16308  and  42800  will 
give  you  the  latter  fix  instead  of  the  former.  Proceed  as  follows 
(assume  the  GRI ' s  have  been  loaded  as  in  Sample  Problem  1): 


DISPLAY  CONTENTS 

USER  RESPONSE 

COMMENTS 

NEXT  OPTION  (Help) ? 
LAT  dd.mmss? 
LONG  ddd.mmss? 
ITD  FOR  9940W? 
ITD  FOR  9940Y? 

7 

36.4736  [RTN] 

121.4658  [RTN] 

1630  8   [RTN] 

42800   [RTN] 

Calibrate. 

36° 47 '36" 

121°45,58M 

The  station  baselines  have  now  been  modified  to 
achieve  calibration.  Now,  reenter  16038  and  42800  into 
the  fixing  routine  to  test  the  calibration.  Calibration 
does  not  modify  the  permanent  station  data  stored  in  the 
HP-75C.  The  calibration  remains  effective  until  the 
GRI's  are  reloaded  by  Option  8.  NOTE:  The  affect  of 
calibration  on  the  accuracy  of  fixes  far  removed  from 
the  benchmark  has  not  been  studied.  The  user  should  use 
the  calibration  option  with  caution. 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


NEXT  OPTION  (Help) ? 
ITD  FOR  9940W? 
ITD  FOR  9940Y? 
LAT   =   36° 47 '35" 
LONG  =  121° 46' 59" 


16308  [RTN] 
42800  [RTN] 

[RTN] 

[RTN] 


Compute  a  fix, 


Lat  of  fix. 
Long  of  fix. 


The  small  discrepancy  between  this  fix  and  the 
benchmark  is  due  to  assumptions  in  the  fixing  algorithm. 
If  higher  precision  is  required,  the  ASF  Correction 
factor  may  be  used  (see  Sample  Problem  5).  Note  that  you 
may  need  to  use  the  Alternate  Solution  Option  2  if  the 
fix  obtained  above  is  near  38°52'57"N  and  116°51'11,W. 


DISPLAY  CONTENTS 
NEXT  OPTION  (iielp)  ? 


USER  RESPONSE 


COMMENTS 

Repeat  with  new  option 
or  press  A  to  quit. 
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HP-75C  SAMPLE  PROBLEM  5 

Suppose  you  are  receiving  an  ITD  of  12153.31  microseconds 
from  9960W  and  44451.83  microseconds  from  9960Y.  These  ITD's  would 
tell  you  that  your  location  is  44°15,05"N  and  67°25,23"W.  Entering 
the  ASF  Correction  Tables  [Ref.  5,6]  with  these  readings,  the  ASF 
Correction  for  GRI  9960W  at  44°15'N  and  67°25'W  is  1.5  micro- 
seconds. Similiarly,  the  ASF  Correction  for  GRI  9960Y  at  these 
coordinates  is  2.7  microseconds.  Enter  these  corrections  into  your 
HP-75C  and  recompute  your  position.  (Load  the  GRI * s  9960W  and 
9960Y  using  Option  8). 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


NEXT  OPTION  (Help)? 
ASF  FOR  9960W? 
ASF  FOR  9960Y? 
NEXT  OPTION  (Help) ? 
ITD  FOR  9960W? 
ITD  FOR  9940Y? 
LAT   =   44° 15' 26" 
LONG  =   67° 26  '26" 
NEXT  OPTION  (Help) ? 


1.5 
2.7 

12153 
44451 


9 

[RTN] 

[RTN] 
1 

.31 

.83 

[RTN] 

[RTN] 


[RTN] 
[RTN] 


ASF  Correction, 


Compute  a  fix. 


Repeat  with  new  option 
or  press  A  to  quit. 


These  ASF  Corrections  have  been  entered  into  the 
computer.  The  range  of  their  validity  must  be  checked 
with  the  ASF  Correction  tables  and  reentered  if  they 
change.  NOTE:  If  you  are  using  the  ASF  Correction  Tables 
and  if  you  wish  to  calibrate  to  a  benchmark  position, 
you  should  first  enter  the  ASF  Corrections  for  your 
benchmark  and  then  use  the  calibration  Option  7  (see 
Sample  Problem  4). 
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HP-75C    PROGRAM    LISTING 


18  RE*  LORAN-C  ROUTINE. 
rev  87-88-83.    1336  HOURS 

26  N9=15  0  REN  N9=l  0^  H 
ASTER  STATIONS 
38  OPTION  ANGLE  RADIANS 
0  OPTION  BPSE  1 
40  DEF  FHR(XjH)  =  INT(H* 
X+.5)/N 

50  C8= 1/298. 26  6  A8=-l  * 
Gl$=" 

60  DEC  FNP(T)  =  129.0439 
8/T- . 40758+ . 00064576438* 
T 

70  DEF  FN9('T)  =  2.741297 
9 /T-. 011 482+. 00832774624 
*T 

88  LISP  0  GOTO  126 
98  C*=KEY$  0  IF  C$="  T^ 
EN  90 

108  C*=OPRC$(C$)  0  FOR  C 
=1  TO  LEN(CW)  0  IF  C$=C 
8«C.C1  THEN  D ISP  0  RETU 
RN 

110  NEXT  C  0  GOTO  98 
128  PISP  "NEXT  OPTION  C 
jCHRKMe^elpP'  0  C6* 
=,1234567890W  8  GOSUB  9 
8 

130  ON  C  GOSUB  916.950. 1 
010,1848.1898,998,1286.1 
98,1358,538,148  0  GOTO  1 
28 

146  BISP  "1  -  LORAN  FIX' 
0  BISP  '2  -  ALTERNATE  S 
OLN-  0  DISP  '3  -  INPUT  H 
EST  LPT/LONG1 
156  DISP  '4  -  HEAD  %   BIS 
T  TO  BEST"  0  DISP  '5  -  P 
REPICT  ITDs- 
166  DISP  '6  -  THO  POINT 
HEAD  %   BIST'  0  DISP  "7  - 
CALIBRATE" 

176  DISP  -8  -  NEH  STP.  I 
Ds-  0  BISP  -9  -  ASF  Corr 
ections1  0  DISP  *A  -  STO 
P- 


188  RETURN 

182  REM  OPT  8 

185  IF  G !*<>'■  THEN  RETU 

RN 

198  F8=l  0  N=l  0  INPUT  • 

1st  GPP  Ml 

288  GOSUB  388  0  GOSUB  39 

8  0  IF  F8=2  THEN  198 

210  F8=l  0  N=2  0  INPUT  ' 

2nd  GRP  ",  1$  0  GOSUB  3 

88 

220  IF  fii=!9*  THEN  GOSUB 
478  ELSE  248 

238  IF  F8=l  THEN  258 

248  GOSUB  398  8  IF  F8=2 

THEN  218 

250  Fl=!  0  F2=l  0  IF  Gl* 

=G2$  THEN  388 

268  IF  P(L1)=P(L2)  AND 
L(1,1)=L(1,2)  THEN  316 

276  F2=-l  0  IF  P(1,1>=P< 

2,2)  AND  L(1,1)=L(2,2)  T 

HEN  318 

288  Fl=-1  0  IF  P(2,1)=P( 

2.2)  AND  L(2,1)=L<2,2)  T 

HEN  316 

296  F2=l  0  IF  P(2,D=P(1 

•  2)  AND  L(2,1)=L<I..2)  TH 
EN  310 

380  DISP  'E.  NO  TRIPLET 
-,G1I;"  t  \;G2$  0  GOSUB 
1370  0  GOTO  198 
310  FOP  1=1  TO  2  0  P1=P( 
1,1)  0  L1=L(1..I)  0  P2=P( 
2,1)  0  L2=L(2,I) 
328  DISP  0  GOSUB  338  0  N 
EXT  I  0  G(l)=8  0  G(2>=8 
8  RETURN 

338  GOSUB  658  0  B(I)=S8 
0  Z<1,I>=Z1  0  Z(2,D=Z2 
348  T=21?82.3593*S8 
350  IF  T>=537  THEN  P=FNP 
(T) 

368  IF  T<537  THEN  P=FNO>: 
T) 
376  T(I)=T+P  8  RETURN 
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HP-75C    PROGRAM    LISTING     (cont.) 


380  L=LEN(I$>  0  ft*=I*[l, 

L-l]  0  B*=I$tL,U  0  B=NU 

H(B$>  0  B$=CHR$(B-32*(B> 

90  >>  0  RETURN 

390  N8=l  0  RESTORE 

400  READ  I9$,N0,P6,L«  0 

*%='• 

410  FOR  1=1  TO  N0  0  READ 
W0*,B9(I),Pl(I),LHn  0 
K$=Mtt«0$  0  NEXT  I 

420  IF  I9$fp$  THEN  450 

430  GOSUB  470  0  IF  F0=2 

THEN  468 

440  RETURN 

45&  N8=N8+1  0  IF  N8<=N9 

THEN  400 

460  BISP  It;'  IS  NOT  CPT 

QLOGEP"  0  F0^2  0  GOSUB  1 

370  0  RETURN 

470  FOP  1=1  TO  N0  0  IF  M 

*[I,I]=B*  THEN  496 

480  NEXT  I  0  F0=2  0  RETU 

PN 

490  IF  N=l  THEN  Gl$=!$ 

500  IF  N=2  THEN  G2*=I* 

510  B(N>=B9(I>  0  P(1,H)= 

FW(RflD<FNH<Pe>))  0  L(L 

N>=RRB<FNH(L0>) 

520  P(Z,H)=FHfi(RflD(FNH(P 

i(I))»  0  L(2,N)=RflD(FNH 

(LKI)))  0  RETURN 

530  DISP  0  PUT  CHR*(13) 

0  END 

540  REH  DIR  SOLN 

550  Z8=SIH(Z1)  0  Z9=C0S< 

ZD  0  P8=SIN(P1)  0  P9=C0 

S(P1)  0  H=-Z8*P9  0  C1=C0 

*fl  0  C2=C0*(l-**t1>/4 

560  B=U-C2)*(1-C2-C1**> 

0  P=C2*<l+Cl*«/2)/B  0  N 
=P9*Z9 

570  S1=8NGLE(P8,N>  0  B8= 
S0/B  0  U=?*(SI-B8)  0  H=l 
-2*P*C0S(U)  0  V=COS(U+B0 
) 


580  X=C2*C2*SIH(B0)*COS< 

B«*(2»V*V-!))  0  Y=2*P*V* 

U*SIN(B8> 

590  S2=B0+X-Y  0  S8=SIN<S 

2)  0  S9=C0SfS2>  0  K=S8R< 

H**-KH*S9-P8*S8)A2> 

600  P2=flTN((P8*S9+N*S8)/ 

K) 

610  S3=RNGLE(P9*S9-P8*S8 

*Z9,-S8*Z8> 

620  H=C1*(1-C?)*S2-C!*C2 

*S8*C0S(S1+S1-S2)  0  L2=L 

1+S3-H 

630  Z2=ANGLE'-(N*S9-P8*S 

8)»-H)  0  RETURN 

648  REH  REVERSE  SOLN 

650  L3=L2-L1  0  P3=(P2-P1 

)/2  0  P4=(Pl+P2)/2 

660  P*=SIN(P3>  0  P7=C0S( 

P3'>  0  P8=SIN<P4)  0  P9=CC 

S(P4) 

A70  H=P7*P7-P8*P8  0  L=P6 

*P6+H*SIN(L3/2)A2  0  B8=2 

*flSIN(SeR(D) 

680  U=2»P8*P8*P7»P7/(1-L 

)  %  v=2*P6*P6*P9*P9/L  0 

X=u+V  0  Y=U-V  0  T=D6/SIN 

(JW)  p  p=4*T*T 

690  E=2*COS(B0)  0  A=M 

0  C=T-(fl-E)/2  0  H1=X*<"P+ 

C*X>  0  B=B+D  0  N2=Y*^B+E 

*Y>  0  N3=B*X*Y 

700  D2=CT*C«*(H1-H2+N3)/ 

64  0  Bl=C8*(T*X-Y)/4  0  S 

8=(T-B1+B2)*SIN(B0) 

710  H=32*T-(2S*T-R)«X-<B 

+4)*Y 

720  F=Y+Y-E*(4-X)  0  G=C8 

*(T/2+C0*«/64)  0  Q=-F*G* 

TBN(L3>/4 

738  L4=(L3^)/2  0  L8=SIN 

(L4)  0  L9=C0S(L4) 

740  T1=ANGLE(P9*L8,P6*L9 

)  0  T2=RNGLE(P8*L8,-P7*L 

9) 
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HP-75C    PROGRAM    LISTING     (cont.) 


758  T9=PI+PI  8  Z1=H0B(T1 

♦T2,T9>  8  Z2=H0D(T1-T2,T 

9)  8  RETURN 

768  REM  FIXING  ROUTINE 

77*  Al=Fl*SIN(fl(l>>  8  Bl 

=C0S<fl(l))-C0S(B(D)  8  C 

1=SIN(B(D) 

788  R2=F2*SIN<fl<2>>  8  B2 

=C0S(P(2))-C0S(B(2))  8  C 

2=SIN(B(2>) 

798  E1=Z<1,1>  8  IF  Fl=-i 
TNEN  E1=Z(2..1> 

888  E2=Z(1,2)  8  IF  F2=-i 

THEN  E2=Z(2..2> 
818  C=P!*C2*C0S(E21-B?*C 

l*COS(El>  8  S=B1*C2*SIN( 

E2>-B2*C!*SIN<E1) 

828  K=82*fll-81*A2  8  P=SQ 

R(OOS*S>  8  G=PNGLE(CS 

) 

838  Z=G+fl8*RC0S<K/R>  8  S 

8=flNGLE<C2*C0S(Z-E2>+fl2. 

B2> 

848  IF  F2=l  THEN  P1=P<L 

2)  8  L1=LU,2> 

858  IF  F2=-l  THEN  P1=P(2 

,2)  8  L!=L<2,2) 

868  Z1=Z  8  GOSUB  558  8  P 

8=P2  8  L8=L2  8  P=flTNaflN 

(P8)/(!-C8)) 

878  P=BEG(P)  8  L=H0B(BEG 

(L8>,368>  8  IF  L>188  THE 

N  L=L-368 

888  P*=FND*<P>  8  L$=FNF* 

(L) 

898  BISP  -LRT  =  ";P$  8  G 

OSUB  1378  8  DISP  "LONG  = 

Ml  9  GOSUB  1378  8  RET 
URN 

988  REM  OPT  1  -  COMPUTE 
FIX 

918  GOSUB  185  8  FOP  1=1 
TO  2  8  GMG1$  8  IF  1=2  T 
HEN  G*=G2* 


928  DISP  -ITD  FOR  -;G$; 
8  INPUT  ■?  -.;  p  8  fl=fl+G< 
!)-B(I>-T(I)  8  IF  RBS(R) 
<T(I)  THEN  948 
938  BISP  *E.  ITD  NOT  VPL 
ID  FOR  *;G*  8  GOTO  928 
948  H(I>=8/21295.8736  8 
NEXT  I  8  GOSUB  778  8  PET 
URN 

958  IF  08=1  THEN  08=- 1  8 
GOTO  965 
%8  08=1 

965  IF  Gl$="  THEN  918 
978  GOSUB  778  8  RETURN 
988  REN  OPT  6 
998  INPUT  -ORIGIN  LOT  dd 
.ftftss?  V;  P  8  INPUT  "OR I 
GIN  LONG  ddd.wss?  ";L 
1888  P8=FNM(PflD(FNH(Pm 
8  L8=RflD'rFNH(L))  8  GOSU 
B  1818  8  GOSUB  1848  8  RE 
TURN 

1818  INPUT  -BEST  LRT  66. 

w»ss?  -•  P  8  INPUT  "BEST 

LONG  ddd.fwss^  ";L 

1828  P5=FNP(R0D(FNHfP>^ 

8  L5=RflD(FNH(D)  8  RETU 

PN 

1838  REH  OPT  4 

1848  P1=P8  8  L1=L8  8  P2= 

P5  8  L2=L5  8  GOSUB  658 

1858  Z=ttOB(BEG(Zl),368> 

8  Z*=FND$(Z) 

1868  DISP  "HEADING  =  ";Z 

t  8  G0SU8  1378  8  DISP  'H 

ISTPNCE  =  -;FNR(3443.917 

♦S8.188);-  n.«i.- 

1878  GOSUB  1378  8  RETURN 

1888  REH  OPT  5 

1898  GOSUB  185  8  GOSUB  1 

178  8  FOR  K=l  TO  2  8  19= 

KKHDOf) 

1188  P2=P(2,K>  8  L2=L(2. 

K)  8  GOSUB  658  8  1=3  8  G 

OSUB  348  8  I9=I9+T(3> 
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HP-75C    PROGRAM    LISTING     (cont.) 


1118  P2=PU,IO  0  L2=L<1, 
K>  0  GOSUB  650  0  1=3  0  G 
OSUB  340  0  I9=I9-T(3) 
1120  I$=G1$  0  IF  K=2  THE 
N  I$=G2$ 

1130  MSP  TIP  FOR  -;lti 
'  =-;FW?a9,100),CHRf<12 
).:-S-  0  COSUB  1378  0  HEX 
T  If 

1148  DISP  CHR*<206);-P¥ 
Lat/Lon  or  V;CHRf(218>:- 
eturn  to  wnu*>"  0  M«='H 
R"  9  GOSUB  98 
1158  IF  C=l  THFH  1890 
1168  RETURN 

1178  INPUT  -LPT  dd.mss? 
■;  PI  0  INPUT  "LONG  ddd 
.wiss?  Ml 

1180  P1=FHP(RR1KFNH(P1>> 
)  0  L1=#MKFW(L1))  0  RE 
TUPH 

1198  RE*  OPT  7 
1208  GOSUB  185  0  GOSUB  1 
170  0  FOR  K=l   TO  2  0  I*= 

Git  0  IF  K=2  THEH  I$=G2* 

1218  DISP  -ITU  FOR  M$; 

0  INPUT  ■?  ■;    19  0  19=1 

9-DCP 

1220  P2=P(2.K>  0  L2=L(2. 

K)  0  GOSUB  650  0  1=3  0  G 

OSUB  340  0  !9=I9-T(3) 

1238  P2=P(1,K>  0  L2=L(1.. 

K>  0  GOSUB  650  0  1=3  0  G 

OSUB  340  0  T(K>=I9+T(3)+ 

G(K) 

1240  NEXT  r  0  RETURN 

1250  BEF  FNH(X) 

1260  S=SGN(X)  0  X=RBSW 
0  H=IMT(X)  0  V=FP(X)*1B 

8  0  FHH=S*<<180*FP(V)/60 

♦INT<V)>/6«+H> 

1270  END  DEF 

1288  DEF  FNfl(X)  =  ATN((1 

-C0)*TRN(X)) 

1290  DEF  FNM(X) 


1388  C$='  •  0  IF  X<8  THE 
H  C*='-- 

1310  X=flBS(X>+ 1/7208  0  X 

8=INTfX)  0  Ci=C$tSTR$<X0 

)iCHR$(l) 

1320  X=*0*(X-X0>  0  X0=IN 

T(X)  0  X*=STR*<100+X8^  0 

cf=mx$r2,3H"- 

1330  X=60*(X-X0>  0  X8=IN 
T(X)  0  X*=STWUW*X0>  0 

FNW=C»X*[2.3ttCHPt(34 
) 

1348  END  DEF 

1358  GOSUB  185  0  DISP  -G 
SF  FOR  ■;Gl$;  0  INPUT  ■? 

•;G(1) 
1368  DISP  'QSF  FOR   -;G2$ 
;   0   INPUT    '?   ■;    G(2V'   0  P 
ETURH 

1378  C*=KEY*  0  IF  C$=" 
THEN  !3?0 
1380  DISP  0  RETURN 
5880  WTfl  4998,2,16.4443 
95, 169.38312, "X", 11808.2 
0.144916,155.53097 
5018  &QTP  -V-,29008.28.2 
34177,178.17302 
5020  DPTP  5939,2,46,4827 
199,67.5537713. "X-, 11080 
,41.151193,69.583909 
5830  WTR  -Y-,25880,46.4 
63218,53.102816 
5040  BflTfi  5970,3,36.1105 
797, -129. 2827279, -HV.  110 
00,42.4437104.-143.43092 
45 

5050  BRTP  -X-, 31008, 35.8 
223871. -126. 322674 1,'Z', 
420*0,26.3624975,-128.08 
56445 

5060  WTP  5990,3,51.5758 
78, 122.220224, -X", 11000, 
55.2620851, 131. 1519648 
5070  WTO  -Y*, 27000, 47.0 
34799. 119. 443953, -Z-, 418 
00,50.3629731,127.212904 
3 
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HP-75C    PROGRAM    LISTING     (cont.) 


5888  BPTP  7938,3,59.5917 

27,45. 1827*7.  "HM 1888, 6 

4.542658,23.552175 

5898  BPTP  "X". 21888, 62.! 

75964, 7.8426538, "Z", 4389 

8,46.463218,53.182816 

5188  BPTP  •♦7938", 3, 24.1 

787888,-153.5853232, -XV. 

11888,42.4437184,-143.43 

89245 

5118  BPTP  -Y", 38888, 26. 3 

62*975,-128.8856445, -Z'. 

49888,9.3245789,-138.895 

497 

5128  BPTP  7968,2,63.1942 

814, 142.48319. -X-, 11888- 

57.262821,152.2211225 

5138  BPTR  -Y-,26888,55.2 

628851,131.1519648 

5148  BPTR  7978,4,62.1759 

64.7.8426538, -H-, 26888, 5 

4.4829872,-8.1736312 

5158  BPTP  -X", 11888,68.3 

88615, -14. 2747. "Y-, 46888 

,64.5426.58,23.552175 

5168  BPTR  "Z", 68888, 78. 5 

45261.8.435869 

5178  BPTR  7988-4,38.5938 

74,85. 1889385- "H-, 11888, 

38.4333818,98.49436 

5188  BPTR  "X", 23888, 26. 3 

155886,97.5888893, -Y", 43 

888,27.8158393,88.865342 

9 

5198  BPTR  "Z", 59888, 34. 8 

346881,77.5446654 

5288  BPTR  7998,3,38.5228 

587, -16. 4386159, -X", 1188 

8,35.3128787,-12.3138245 

5218  BPTR  -Y-,29888,48.5 

8?895, -27. 528152, *Z' ,478 

88,42.8336515,-3.1215512 


5228  BPTR  8978-3,39.5187 
54, 87. 291214, -H", 11888, 3 
8.593874,85.1889385 
5238  BPTP  -X", 28888. 42. 4 
258683, 76. 4933862, -Y-, 44 
888,48.3649844,94.331846 
9 

5248  BPTP  9948,3,39.3386 

621,118.495637 

5258  BPTP  •«*, 11888,47.8 

34799, 119.443953, -XV.  278 

88,38.465699,122.2944529 

5268  WTP  -YV. 48888, 35.1 

91818,114.4817435 

5278  BPTR  9968,4,42.4258 

683, 76. 4933862, -H", 11888 

,46.4827199,67.5537713 

5288  BPTP  -X*, 25888, 41.1 

51193,69.583989, "Y-, 3988 

8,34.8346881,77.5446654 

5298  WTP  'Z", 54888, 39. 5 

18754-87.291214 

5388  BPTP  9978,4,24.4883 

597, -141. 1938383, -W,  118 

88,24.1787888,-153.58532 

32 

5318  BPTR  "X", 38888, 42. 4 

437184.-143.4389245. 'YV. 

55888,26.3624975,-128.88 

56445 

5328  BPTR  "Z", 75888, 9. 32 

45789.-133.895497 

5338  BPTP  9998,3,57.8912 

265, 178.1 586789, -X*. 1188 

8,52.494484,-173.1848974 

5348  BPTP  -Y-, 29888, 65.1 

448386, 166.531255, -Z-, 43 

888.57.262821,152.221122 

5 
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III.  THE  SHARP  PC-1500  (RADIO  SHACK  TRS- 80  PC-2) . 

PC-1500  (TRS-80  PC-2)  USER  INSTRUCTIONS. 

The  options  are  labeled  1  through  9,  A  or  H.  The  H  (Help) 
option  will  review  options  1  -  9  &  A.  To  step  through  the  HELP 
menu,  use  the  ENTER  key.  Option  6  (TWO  POINT  HEAD  &  DIST)  is  a 
stand-alone  option  which  may  be  used  by  itself.  For  LORAN  fixing, 
Option  8  must  be  selected  first  to  load  the  station  GRI's.  If 
Option  8  is  not  used,  then  any  option  requiring  the  GRI's  will 
automatically  enter  Option  8  first.  Option  1  (LORAN  FIX)  must  be 
selected  at  least  once  before  Option  2  (ALTERNATE  SOLN)  can  be 
selected.  Once  the  GRI's  are  loaded,  any  of  the  remaining  options 
may  be  used  at  any  time.  When  first  run,  the  ASF  Correction 
Factors  are  set  to  zero;  if  they  are  subsequently  entered  as 
non-zero  values  they  will  remain  set  until  either  changed 
manually,  by  selection  of  Option  9,  or  until  the  program  is  RUN 
once  more.  In  all  cases,  South  latitudes  and  East  longitudes  must 
be  negative  for  input  and  will  be  labeled  by  a  minus  sign  '-'  on 
output.  For  ease  of  reference  the  options  are  listed  by  number 
rather  than  by  sequence  of  selection. 


INSTRUCTION 


DISPLAY 


INPUT 


Run  Program 


RUN. 10 


LORAN  FIX 
ITD  for  1st  GRI 
ITD  for  2nd  GRI 
Lat  of  fix 
Long  of  fix 


NEXT  OPTION  (Help) ? 

ITD  FOR  gri  1? 
ITD  FOR  gri  2? 
LAT  =  dd°mm'ss" 
LONG  =  ddd°mm'ss" 
NEXT  OPTION  (Help) ? 


1st  ITD  [ENTER] 
2nd  ITD  [ENTER] 

[ENTER] 

[ENTER] 


2.  ALTERNATE  SOLN 
Lat  of  fix 
Long  of  fix 


LAT  =  dd°mm'ssH 
LONG  =  ddd°mmlss" 
NEXT  OPTION  (Help)? 


[ENTER] 
[ENTER] 


3.  INPUT  DEST  LAT/LONG 


DEST  LAT  dd.mmss? 
DEST  LONG  ddd.mmss? 
NEXT  OPTION  (Help) ? 


Lat  of  dest. 

[ENTER] 
Long  of  dest 

[ENTER] 


4.  HEAD  &  DIST  TO  DEST 

Heading  to  dest.     HEADING  =  dd°mm'ss" 
Dist.  to  dest.      DISTANCE  =  mmm.ff  n.mi 

NEXT  OPTION  (Help) ? 


[ENTER] 
[ENTER] 
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PC-1500  (TRS-80  PC-2)  USER  INSTRUCTIONS  (cont.) 


INSTRUCTION 
5.  PREDICT  ITDs 

1st  ITD 
2nd  ITD 


6.  TWO  POINT  HEAD  & 


Heading  to  dest. 
Dist.  to  dest. 


7.  CALIBRATE 


DISPLAY 

LAT  dd.mmss? 
LONG  ddd.mmss? 
ITD  FOR  gri  1  = 

xxxxx.xx  ys 
ITD  FOR  gri  2  = 

xxxxx.xx  ys 
Hew  Lat/Lon  or 
Return  to  menu 

NEXT  OPTION  (Help) ? 

DIST 

ORIGIN  LAT  dd.mmss? 

ORIGIN  LONG  ddd.mmss? 

DEST  LAT  dd.mmss? 

DEST  LONG  ddd.mmss? 

HEADING  =  dd°mm,ss" 
DISTANCE  =  mmm.ff  n.mi 
NEXT  OPTION  (Help) ? 


LAT  dd.mmss? 
LONG  ddd.mmss? 
ITD  FOR  gri  1? 
ITD  FOR  gri  2? 
NEXT  OPTION  (Help) ? 


8.  NEW  STA.  IDs 


9.  ASF  Corrections 


1st  GRI? 
2nd  GRI? 
NEXT  OPTION 


(Help)? 


ASF  FOR  gri  1? 
ASF  FOR  gri  2? 
NEXT  OPTION  (Help) ? 


INPUT 


Lat   [ENTER] 
Long  [ENTER] 

[ENTER] 

[ENTER] 

N  repeats  Option  5 
R  goes  to  next  line 


Lat  of  origin 

[ENTER] 
Long  of  origin 

[ENTER] 
Lat  of  dest. 

[ENTER] 
Long  of  dest. 

[ENTER] 

[ENTER] 

[ENTER] 


7 
Lat  Of  fix   [ENTER] 
Long  of  fix  [ENTER] 
1st  ITD  [ENTER] 
2nd  ITD  [ENTER] 


8 


Sta 
Sta 


1st 
2nd 


ID 
ID 


ASF 
ASF 


[ENTER] 
[ENTER] 


[ENTER] 
[ENTER] 


A.  STOP 
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PC-1500  (TRS-80  PC-2)  SAMPLE  PROBLEM  1 


You  are  sailing  off  of  the  coast  of  California.  Your 
destination  is  Moss  Landing  at  about  36°48'N  and  121°47'W.  Your 
present  ITD  from  station  pair  GRI  9940W  is  16019  and  your  ITD  from 
station  pair  GRI  9940Y  is  42585.  What  is  your  current  position  and 
what  is  the  heading  and  distance  from  your  current  position  to 
Moss  Landing?  We  assume  that  the  GRI ' s  9940W  and  9940Y  have  not 
yet  been  loaded. 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


NEXT  OPTION  (Help) ?  3 

DEST  LAT  dd.mmss?         36.48  [ENTER] 

DEST  LONG  ddd.mmss?      121.47  [ENTER] 


Input  destination. 
36°48' 
121°47  ' 


NEXT  OPTION 
1st  GRI? 
2nd  GRI? 


(Help)? 


9940W 
9940Y 


8 

[ENTER] 

[ENTER] 


Input  station  ID's, 


NEXT  OPTION  (Help) ? 
ITD  FOR  9940W? 
ITD  FOR  9940Y? 
LAT  =  39  14  19 
LONG  =  115  50  52 


1 

16019 

42585 

[ENTER] 

[ENTER] 


Compute  a  fix. 


Latitude  of  fix. 
Longitude  of  fix, 


This  fix  is  in  Nevada.  There  are  always  two 
solutions  to  the  LORAN  equations.  The  navigator 
must  select  the  correct  one. 


NEXT  OPTION  (Help)? 
LAT  =  35  00  01 
LONG  =  125  00  09 


[ENTER] 
[ENTER] 


Select  alternate  soln, 
Latitude  of  fix. 
Longitude  of  fix. 


This  fix  is  in  the  Pacific.  The  navigator 
knows  it  is  the  correct  fix  because  of  his 
estimate  of  where  he  should  be. 


NEXT  OPTION  (Help)?  4 

HEADING  =  54  34  11  [ENTER] 

DISTANCE  =  190.38  n.mi.     [ENTER] 


Find  heading  &  distance 
to  destination. 
Heading  to  Moss  Landing. 
Distance  to  Moss  Landing 


NEXT  OPTION  (Help) ? 


Repeat  with  new  option 
or  press  A  to  quit. 
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PC-1500  (TRS-80  PC-2)  SAMPLE  PROBLEM  2 


How  far,  and  in  what  direction,  is  Corvallis,  Oregon 
(44°34*N,  123°16'W)  from  Cupertino,  California  (37°19'N, 
122°02'W)? 


DISPLAY  CONTENTS 

NEXT  OPTION  (Help) ? 

ORIGIN  LAT  dd.mmss? 
ORIGIN  LONG  ddd.mmss? 
DEST  LAT  dd.mmss? 
DEST  LONG  ddd.mmss? 
HEADING  =  353  02  59 
DISTANCE  =  43  8.32  n.mi 


USER  RESPONSE 

COMMENTS 

6 

Two 

point  distance 

and 

heading. 

37.19  [ENTER] 

37°19' 

122.02  [ENTER] 

122°02' 

44.34  [ENTER] 

440341 

123.16  [ENTER] 

123°16' 

[ENTER] 

[ENTER] 

NEXT  OPTION  (Help) ? 


Repeat  with  new  option 
or  press  A  to  quit. 


PC-1500  (TRS-80  PC-2)  SAMPLE  PROBLEM  3 


Suppose  that  you  want  to  know  what  ITD's  you  would  expect  to 
receive  at  35°  North,  125°  West  and  at  36°27'  North,  126°54'  West 
from  9940W  and  9940Y.  To  determine  these  ITD's,  proceed  as  follows 
(assume  that  the  GRI ' s  have  been  loaded  as  in  Sample  Problem  1): 


DISPLAY  CONTENTS 

NEXT  OPTION  (Help) ? 

LAT  dd.mmss? 

LONG   ddd.mmss? 

ITD   FOR   9940W    =    16019.35 

ITD   FOR   9940Y   =    42584.71 

Hew  Lat/Lon  or 

Return  to  menu? 
LAT  dd.mmss? 
LONG  ddd.mmss? 
ITD  FOR  9940W  =  15572.32 
ITD  FOR  9940Y  =  43006.15 
Hew  Lat/Lon  or 

Eeturn  to  menu? 


SER  RESPONSE 

COMMENTS 

5 

35  [ENTER] 

125  [ENTER] 

[ENTER] 

[ENTER] 

Predict  ITD's 
35° 
125° 

N 

36.27  [ENTER] 

126.54  [ENTER] 

[ENTER] 

[ENTER] 

36° 27  ' 
126° 54' 

R 


NEXT  OPTION  (Help) ? 


Repeat  with  new  option 
or  press  A  to  quit. 
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PC-1500  (TRS-80  PC-2)  SAMPLE  PROBLEM  4 


Suppose  you  are  receiving  an  ITD  of  16308  from  9940W  and 
42  800  from  9940Y.  These  ITD's  would  tell  you  that  your  location  is 
36°47,55"N  and  121°47'llnW.  However,  you  know  that  your  position 
is  benchmarked  to  be  at  36047'36"N  and  121°46'58"W,  and  you  wish 
to  calibrate  your  PC-1500  so  that  the  ITD's  of  16308  and  42800 
will  give  you  the  latter  fix  instead  of  the  former.  Proceed  as 
follows  (assume  the  GRI's  have  been  loaded  as  in  Sample  Problem 
1): 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


NEXT  OPTION  (Help) ? 
LAT  dd.mmss? 
LONG  ddd.mmss? 
ITD  FOR  9940W? 
ITD  FOR  9940Y? 


36.4736  [ENTER] 

121.4658  [ENTER] 

16308  [ENTER] 

42  800  [ENTER] 


Calibrate. 

36047'36" 

121°45,58n 


The  sta 

ach 

ieve  cali 

the 

fixing  r 

does  not  mod 

pc-: 

1500.  The 

GRI 

's  are  re 

cal 

ibration 

the 

benchmar 

the 

calibrat 

tion  baselines  have  now  been  modified  to 
bration.  Now,  reenter  16038  and  42800  into 
outine  to  test  the  calibration.  Calibration 
ify  the  permanent  station  data  stored  in  the 

calibration  remains  effective  until  the 
loaded  by  Option  8.  NOTE:  The  affect  of 
on  the  accuracy  of  fixes  far  removed  from 
k  has  not  been  studied.  The  user  should  use 
ion  option  with  caution. 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


NEXT  OPTION  (Help) ? 
ITD  FOR  9940W? 
ITD  FOR  9940Y? 
LAT   =   36  47  35 
LONG  =  121  46  59 


16308  [ENTER] 

42  800  [ENTER] 

[ENTER] 

[ENTER] 


Compute  a  fix 


Lat  of  fix. 
Long  of  fix. 


The  small  discrepancy  between  this  fix  and  the 
benchmark  is  due  to  assumptions  in  the  fixing  algorithm. 
If  higher  precision  is  required,  the  ASF  Correction 
factor  may  be  used  (see  Sample  Problem  5).  Note  that  you 
may  need  to  use  the  Alternate  Solution  Option  2  if  the 
fix  obtained  above  is  near  38052'57KN  and  116°51,11,W. 


DISPLAY  CONTENTS 
NEXT  OPTION  (Kelp)? 


USER  RESPONSE 


COMMENTS 

Repeat  with  new  option 
or  press  A  to  quit. 
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PC-1500  (TRS-80  PC-2)  SAMPLE  PROBLEM  5 


Suppose  you  are  receiving  an  ITD  of  12153.31  microseconds 
from  9960W  and  44451.83  microseconds  from  9960Y.  These  ITD's  would 
tell  you  that  your  location  is  44°15'05MN  and  67025'23nW.  Entering 
the  ASF  Correction  Tables  [Ref.  5,6]  with  these  readings,  the  ASF 
Correction  for  GRI  9960W  at  44°15'N  and  67°25'W  is  1.5  micro- 
seconds. Similiarly,  the  ASF  Correction  for  GRI  9960Y  at  these 
coordinates  is  2.7  microseconds.  Enter  these  corrections  into  your 
PC-1500  and  recompute  your  position.  (Load  the  GRI • s  9960W  and 
9960Y  using  Option  8)  . 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


NEXT  OPTION  (Help) ? 
ASF  FOR  9960W? 
ASF  FOR  9960Y? 
NEXT  OPTION  (Help)? 
ITD  FOR  9960W? 
ITD  FOR  9940Y? 
LAT   =   44  15  26 
LONG  =   67  26  26 
NEXT  OPTION  (Help) ? 


1.5  [ENTER] 

2.7  [ENTER] 

1 

12153.31  [ENTER] 

44451.83  [ENTER] 

[ENTER] 

[ENTER] 


ASF  Correction, 


Compute  a  fix. 


Repeat  with  new  option 
or  press  A  to  quit. 


These  ASF  Corrections  have  been  entered  into  the 
computer.  The  range  of  their  validity  must  be  checked 
with  the  ASF  Correction  tables  and  reentered  if  they 
change.  NOTE:  If  you  are  using  the  ASF  Correction  Tables 
and  if  you  wish  to  calibrate  to  a  benchmark  position, 
you  should  first  enter  the  ASF  Corrections  for  your 
benchmark  and  then  use  the  calibration  Option  7  (see 
Sample  Problem  4) . 
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PC-1500  (TRS-80  PC-2)  PROGRAM  LISTING 


10: REM  LORAN-C  R 
OUTJNE.  ^ev  06 
-23-83.  1730  H 
OURS. 

20:N9=15:REf1  N9= 
n  OF  MASTER  ST 
ATJONS 

25:D1M  W$<4)*1,  D9 
(4),  PI  (4),  LI  (4 
),  G$C2)*6,  D<2) 
,  PC2,  2),  L(2,  2) 
,  BC2),  TC3),  ?(2 
,  2 ) ,  A  ( 2 ) 

26: DIM  GC2) 

30:RAD]AN  :C0=l/2 
98.26: A0=~1:RD 
=  P1    /180:G1$=' 
::GOTO  120 

40:X=1NT  CM*X+.5) 
/M: RETURN 

50:P=129.04398/T- 
. ^0758+. 645764 
38E-3*T: RETURN 

60:p=2. 7412979/T- 
.01 1402+. 32774 
624E-3*T 
C$=INKEY$  : IF 
C$='"'THEN  70 


140; 


70 
80:  C 


90 


100 
120 


130 


CHR$  <C-32*(C> 
90)) 

FOR  C=1T0  LEN 
(C0$): IF  C$  = 
MID$  CC0$,  C,  1) 
THEN  RETURN 
NEXT  C:GOTO  70 
PAUSE  "NEXT  OP 
TION  (He lp)?": 
C0$=" 123456789 
AH":GOSUB  70: 
PAUSE  "  " 
ON  CGOSUB  910, 
950,  1010,  1040, 
1090,  990,  1200, 
190,  1320,  530,  1 
40:GOTO  120 


150; 


160: 


170 


PRINT  ' 1 

-  LOR 

AN  FIX'  :  i 

DR1NT 

"2  -  ALTERNATE 

SOLN":pF 

*INT  " 

3  -INPUT 

DEST 

LAT/LONG 

■ 

PRINT  "4 

-  HEA 

D  &  DIST 

TO  FI 

X|::PRINT 

'  5  - 

PREDICT 

I  TDs  ,; 

180 
185 

190: 


200 
210 

220 

240 
250 
260 

270 


PRINT  '6  -  TUO 
POINT  DIST  & 

HEAD':PRINT  "7 
CALIBRATE1 

DR1NT  '8  -  NEui 
STA.  IDs'. 

PR1NT  "9  -  ASr 

Correct  i  ons  '' : 

PRINT  ,:A  -  STO 
p.. 

RETURN 

IF  Gl$<  >'  '  THEN 
RETURN 

F0=j :M=l; INPUT 
'1st  GRI°  ' ; 1$ 
:  PAUSE  ,; 
COSUB  380: 
COSUB  390; 
0-2THEN  130 
P0-1 :N^2: INPUT 
"2nd  GR1?  '  ;  1$ 
: PAUSE  '    : 
GOSUB  380 
IF  A$=19$THEN 
GOSUB  470: IF  F 
8=1  THEN  250 
GOSUB  390: IF  F 
0-2THEN  210 
Fl  =  l:F2--l  :  IF  G 
1$=62$THEN  300 
IF  PCI,  1 >  =P < 1 , 
2) AND  L<1,  1)=L 
(1, 2)THEN  310 
F2=-l:  IF  PCI,  1 
)=P<2,  2) AND  LC 
1,  1  )=LC2, 2) 
THEN  310 
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PC-1500  (TRS-80  PC-2)  PROGRAM  LISTING  (cont.) 


28e 

290 

300 
310 

320 

330 

340 
350 
360 
370 
380 


390 
400 

410 


420 
430 


Fi=-l:  IF  P(2,  1       440 

)=p(2,  2^AND  L(       450 

2,  1  >=L(2, 2) 

THEN  310  460 

F2=l:  IF  P(2,  1 ) 

=PC1, 2>AND  LC2 

,  1)=L(1,  2) THEN       470 

310 

PRINT  "ERROR: 

NO  TRIPLET.  ■ :       480 

GOTO  190 

FOR  l  =  ]TO  2:pl       490 

=P(1,  I ):L1=L<1 

,  J  );D2  =  D(2,  I  >: 

L2-L  :2, : '  500 

PAUSL  '  L*JC"-k  i  na 

': GOSUB  330: 

NEXT  1:G(1)=0:       510; 

G(2>=0: RETURN 

COSUB  650:  BH  ) 

=S0:Z<1, I >=Z1 • 

7(2*  I )=Z2 

'-21282. 3593*5 

0 

IF  T>=--537THEN 

GOSUB  50 

IF  T<537THEN 

GOSUB  60 

T(]  )  =  TfP: 

RETURN 

L-LEN  (1$^:A$- 

LEFT$  (]$, L-l  ) 

:B$=RI6HT$  <]$ 

,  1 ):B=ASC  (B$) 

:B$=CHR$  (B-32 

*(B>90)):  560 

RETURN 

:N8=1 : RESTORE 
:READ  19$, N0,  p0 

, L0  570 

;FOR  1  =  1  TO  N0: 

READ  U$(I ),  D9( 

I ) ,  p  i  ( i  y ,  l  i  <  i ) 

:NEXT  1 
:  IF  19$<>A$THEN 

450 
: GOSUB  470: IF  F 

0=2THEN  460 


:RETURN 

. N8-N8-1 : Ir  N8< 
-N9THEN  400 

: PAUSE  ]*;   IS 
NOT  CATALOGED' 
;F0^2: RETURN 

;POR  1  =  1  TO  N0: 
IF  UI$<1  )=B$ 
THEN  490 
NEXT  1:^0=2: 
RETURN 

IF  N=1THEN  LET 
Gl$=l$: G$( 1 )-] 
$ 

Ir  N=2"THEN  LE" 
G2$-l$: G$<2)=1 
$ 

D(N)=D9( I ): *=R 
D*DEG  (pp): 
COSUB  1270: p:i 
, N)=X: L( J. N)=R 
D^DEG  (l0) 
520:X=RD*DEG  (PI (J 
>): GOSUB  1270; 
D(2, N)=X: L<2,  N 
>--RD*DEG  (LI  (  I 
) ) : RETURN 
DAUS£  :FND 
REM   DIR  SOLN 


0 
550 


78-S1N  (Zl ^:?S 
=COS  CZ1):°8= 
SIN  (P!):P9= 
COS  (PI ^:n--78 
*P9:C1=C0*M:C2 

=  C0*( i-n*n>/4 

D=(l-C2)*<l-C2 

-Cl*n):P=C2*(l 
+Cl*M/2)/D:N=P 
9*Z9 

YY=N:XX=P8: 
COSUB  1260:S1= 
AN: D0=S0/D:U-2 
*(Sl-D0):W-]-2 
*P*COS  (U):U= 
COS  (U+D0) 
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PC-1500  (TRS-80  PC-2)  PROGRAM  LISTING  (cont.) 


580:X-C2*C2*S1N  CD 
0)*COS  CD0*<2* 
U*U-1 )): v=2*p* 
U*U*SIN  CD0) 

590:S2=D0+X-Y:S8= 
SIN  CS2):S9- 
COS  CS2):K-SQR 

<n*n+f n*S9-P8* 

S8)^2) 

600:P2-ATN  CCP8*S9 
+N*S8)/K) 

610: vy=-S8*Z8:XX-D 
9*S9-P8*S8*Z9: 
GOSUB  1260: S3- 
AN 

620:H-C1*( 1-C2)*S2 
-C1*C2*S8*C0S 
CS1+S1-S2):L2= 
L1+S3-H 

630:  vy=-(1:  XX--<:N*S 
9-P8*S8):G0SUB 
1260:Z2-AN. 
RETURN 

640: REM   REU  SOlN 

650:L3=L2-L1 :d3-CP 
2-Dl )/2: P4=(PJ 
+P2>/2 

660:p6-SlN  CP3):P? 
-COS  CP3):p8- 
SIN  (P4):p3- 
COS  (P4) 

670:H=^P7*P7-P8*P3: 
L=P6*P6+H*SIN 
CL3/2)^2:D0-2* 
ASN  (SQR  CD) 

680:U=2*P8*P8*P7*P 
7/Cl-L):0=2*P6 
*P6*P9*P9/L:X- 
U+U: v=u-U: T-D0 
/SIN  CD0):D-4* 
T*T 

690:E-2*COS  CD0):A 
=D*E:C=T-(A-E) 
/2:N1=X*(A+C*X 
):B=-D+D:N2-Y*C 
B-*-E*Y>:N3-D*X* 
Y 


710 
"20 


700:D2=C0*C0*<N]-N 
2+N3)'64:DJ=C0 
*<T*X-Y)/4:S0- 
CT-D1+D2)*SIN 
(D0) 

M-32*T-C20*T-A 
)*X-C8+4)*Y 
F  =  Y+v-£*(4-X) : 

G=C0*(T/2+C0*n 
/g4  ) :  Q--F*G* 

TAN  (L3)/4 
L4=<L3+Q)/2: L8 
-SIN  CL4):L9= 
CCS  (L4^ 
VV---P6*L3:  XX-D9 
*L8:C0SUB  1260 
:  T  i  ^QN:  vvy  --p  '* 
i_3;  XX-P8*L8: 
CGSUB  1260: T2- 
AN 

750: n-P]  +DI  :X=T- 
-tT2:C0SUB  1250 
: 7]^x; X=T1-T2: 
CQSJB  1250:72= 
X:RETURN 

760:REH   FIXING  RC 
UTINE 
A3-F1*SIN  (AC 


'30 


740 


770 


780 


> >-COS  (B( 1 )) 
CI -SIN  (BCD) 


A2=F2*SIN  (AC 2 
)):B2-C0S  (AC 2 
))-COS  CBC2)): 
C2-SIN  (BC2)) 

790:E1-ZC1,  1): IF  F 
1--1THEN  LET  E 
1=Z<2,  1 ) 

800:E2-ZC1, 2): IF  F 
2=-lTHEN  LET  E 
2-ZC2, 2) 

810:C-B1*C2*COS  (E 
2)-B2*Cl*C0S  ( 
E1):S-B1*C2* 
SIN  (E2)-B2*C1 
*S1N  CE1 ) 
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PC-1500  (TRS-80  PC-2)  PROGRAM  LISTING  (cont.) 


820 


830 


840. 


850 


k=B2*Al-Bl*A2: 

R=SQR  CC*C+S*S 

) : YY=S:XX=C: 

COSUB  1260: G-A 

N 

7=G+A0*ACS  (K/ 

R): YY=B2:XX=C2 

*COS  (Z-E2)-A2 

:GOSUB  1260: S0 

=AN 

IF  r2^1THEN 


LET  P]^D( 
L3=L<1, 2) 


-  j 


2  J 


r2^-l THEN 


_t 


L3=LC2,  2) 

860 

:Z1 =2: GOSUB  550 

:P0^D2: L0^L2: z 

--ATN  (TAN  (P0) 

/(1-C0)) 

870 

:p=P/RD:X=L0/RD 

:n-=360:GOSUB  I 

250:L=X: IP  L>] 

80THEN  LET  L=L 

-360 

880; 

;*^D: GOSUB  :280 

:  p$^  $  .  y-|_  . 

COSUB  1280 

890; 

CR]N~  ' _A~ 

;p$:PRJNT  '  LON 

G  =  ' ;C$: 

RETURN 

900: 

REM   OPT  1 

910: 

;GOSUB  185:POR 

1  =  1  TO  2 

920: 

PAUSE  "ITD  FOR 

' ;G$C1 ); : 

INPUT  "?  ";A: 

PAUSE  : PAUSE  " 

;':A=A+G(1)-D< 

I  )-T< I ):  IF  ABS 

(fi)<T(I )THEN  9 

40 

930: 

DR1NT  "ITD  NOT 

UAL  ID  FOR  ";G 

$(1  ):GOTO  920 

340:  a<  ]  ->^a  '21295.  8 

?36:NEXT  ]. 

COSUB  770: 

RETURN 
950: IF  A0=2THEN 

LET  A0=-l:GOTO 

970 
960: A0=1 
970:GOSUB  770: 

RETURN 
980: REH   Dp~  6 
390: INPUT  'ORIGIN 

L  A"  dd  .  nnms  s  ? 

IN  l 0 N G  d  d  d  ■  in  m 
ss°  '  ■ ; L 

1000:X=RD*DEG    Cp" 
:GOSJB    12^0; 
P0=X: L0=RD* 
DEG    ( L  i  • 
COSUB    1010: 
COSUB    1040: 
RETURN 

1010.  INPJT    '  des' 
_AT    dd.rimss0 
;p; INPUT    ' 
D E S """    L 0 NT    d d 
d .  ttis  ?  '■'       i  _ 

•02C:  v-RC^5l:.   - 
: GOSJB  12?0: 
D5=X:L5=RD* 
DEG  (L): 
RETURN 

1030:REM   OP"  4 

1040:P]=P0:L1=L0: 
P2=P5:L2-L5: 
GOSUB  650 

1050:n-360:X-Zl/R 
D:GOSUB  1250 
: GOSUB  1280 

1060:pRlNT  'HEADI 
NG  =  '■;L$:1^ 
100:X-3443.9 
17*S0:  GOSUB 
^0:PR1NT  "Dl 
STANCE  =  ";X 
;   n  .  on  !  .  ' 

1070: RETURN 


-24- 


PC-1500  (TRS-80  PC-2)  PROGRAM  LISTING  (cont.) 


2  080: REM   0pT  5 
1090. GOSUB  185: 
GOSUB  \\7%: 
-OR  k^lTO  2: 
19=T(k)-D(k) 
I100:p2-P(2, K):L2 
=L<2i K): 
GOSUB  650: 1= 
3: GOSUB  340: 
1 3= I9+T<3) 

I : 10:p2  =  Pf  1, k):L2 
-L<1, K): 
GOSUB  650: i  = 
3;COSUt  3*0: 
•■1-  IS-TC3) 

1  120:^100:  v_;  9. 


!0 


GOSUB  4p 
DRINT  '  V 


OR  ' ; G$Ck) ; 

=  '  ;X:  PAUSE. 

r  ":NEXT  k 

i  ;40 

:  ^AUSE  'New  ! 

at  /  i  onq  o^  R 

p i u^n  t  o  me  n 

:C0$  'NR' 

- GOSUB  70 

:  :50; 

:  :r  C- 1 THEN  ; 

090 

i  160 

: RE "URN 

•  120; 

: INPUT  "LA"  d 

d .  Tl  Tl  .  S  S  ?  ' '  ;  p 

.1  ;  INPUT  '  LON 

G  ddd.  nnms  s? 

: ;L1 :pAUSE  " 

1  180: 

;X=RD*DEG  (PI 

>: GOSUB  1270 

:P]=X:L1^RD* 

DEG  CL1): 

RETURN 

1  190: 

REfl   OPT  7 

1200: 

GOSUB  185: 

GOSUB  1170: 

POR  k-lTO  2 

1210: PAUSE  "I  TO  P 
OR  ' ;G$<K); : 
INPUT  "?  ' , I 
9:PAUSL  : 
PAUSE  '  : : 19 
=J9-DCK) 

1220:p2^P(2,  k)-L2 
=LC2, K): 
GOSUB  650:  I - 
3: GOSUB  34  0: 
!9=J9~T(3) 


=LC 1,  k): 

GOSUB  £50.  I  - 

3: GOSUB  340; 

T(k)^]9-T  3N 

rGCk) 

1240 

:NEXT  k: 

RETURN 

1250 

;x=x-n*!N~  <:>: 

/M): RETURN 

1260 

;  AN=ATN  (vv-  '  : 

XX+1E-3*' XX= 

0)))+Pj   X   ■■ 

<0): RETURN 

1270, 

;y-ATN  CC 1-00 

)*TAN  (X; ; : 

RETURN 

1230, 

; C  $- '   •  '  ■   • 

8THnN  ^l.t  C$ 

-  •  ■  -  ■  :  y  -  -  v 

1290: 

; X-X  +  1  /7200:  X 

O^JNT  (X):S$ 

=C$+STR$  CX0 

)  +  "  ' 

1300: 

; X=60*(X-X0>: 

X0=INT  (X):X 

$=STR$  (100^ 

X0>: G$^C$- 

RIGHTS  <X$,  2 

)  +  "  " 

1310: 

:  X-60*(X-X0) : 

X0^1NT  (X):X 

$=STR$  (100+ 

X0):C$=C$+ 

RIGHTS  (X$,  2 

)  +  i:  ':  RETURN 
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PC-1500  (TRS-80  PC-2)  PROGRAM  LISTING  (cont.) 


•32?:C0SUB  185: 

pAUSc  ' ASr  r 
OR  '■  ;  G$<  i  )  j 
INPUT  '  "  '■;  G 
'.  1  ) : PAUSE  : 
PAUSE  : 

1330:PAUSE  'ASr  F 


OR 


)S  '  <l   . 


INPUT  '"> 
'2): PAUSE  : 
PAUSE    ' : 
RETURN 
S000.  DAT;  "4390'  . 


,  i  g  .  ^  a  4  3 

3 1"'  - 1 2 . 
7, 


_K 


037 

■'  10:  DA  rA   -  ■ 

.  178. 17 
220:  DATA  '  53 .-  ! 
?.  d<=  ^B2" 

v.  ' 
33,  53.  5 
33  • 
2:  Al  •  i  ' 

00,  --  .  -63: 

.  53.  102815 
040:  ATA  '  5370  . 

3,  36. : : 05737 
. -123.202727 
9,  ' 'W,  1  1000, 
42.4437J04,  - 
143. 43032^5 

5050: DATA  ' -X'  , 

00, 35.022387 
1, -126.32267 

4,  '  ?'  ,  ^ 2030, 
26. 3624375, - 
128.0856445 

5060:  DATA  ,;5330". 
3,  51  .5^5878, 
122.220224, " 
X'  ,  1  1000,  55. 
2620851, 131 . 
1513648 


5070: DATA  ' v<  ,  270 
03, 47. 034793 

,113.  443a  -, 
7    . 4]00e, 50 

. 3623731,  12" 
212304  3 
5080: DAT3  ■   7930' , 

3,  59.53172"', 

45. 102747,   y 
•  1  1000,  64.  '- 

42658,  23.  55. 
1  75 
5030. SAT  A  ' V ,  210 
00, 62. \?z  - 
■      .  042553'r- 
?'  ,  4  3000,  £g, 

463213, 53. 10 
~  gig 

0;  HATA  '  i'7330' 
•  3,  74  .  '  •' 
8.  -  1 53. 59537 
32,  v  ,  .  )00 
.  4  2 .  4  4  3  ? :  0  .. 
-143.4303 


51 10: 7ATA 


•  300 


00, 26. 362497 
5,-128.08  ' 
■-  .    .  49002 

■  . 32437?  , 

1 38. 0954  37 
:.v. DATA  '  ^360' . 

2,  63.  13  . 

, 142.4831^. 

X:  ,  1  1000,  S."7. 

262021, 152.2 

21 1225 
5130: DATA  "v-,  260 

00, 55.262085 

1, 131. 151364 

8 
5140:DATA  "?370", 

4, 62. 175964, 

7.0426538,  '  U 

",  26000,  54.4 

829872,  -8.  17 

36312 
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PC-1500  (TRS-80  PC-2)  PROGRAM  LISTING  (cont.) 


5150:DATA  ''X'  ,  1  10 
00, 68. 380615 
,  -14.274?,  v 
'  ,  46000,  64.5 
42658, 23.552 
175 

5160:DATA  "Z", 600 
00,  70.545261 
, 8.435869 

5170: DATA  "7980'  . 
4, 30.593874, 
85. 1009305, 
uJ",  1  1000,  30. 
4333018, 90. - 
94  36 

5180: DATA  'X'  ,  230 
00, 26. 315500 

6,  97.500009:- 
.  "Y",  43000,  2 

7.  0158393,  80 
.0653429 

5I90:DATA  "?'  ,  590 
00,  34. 034609 
1 , "7. 5446654 

5200: DATft  '  ^990' , 
3,  38.5220537 
,  -16.4306:5° 
.  X'  .  1  1000.  ? 
5. 3120787, -1 
. 3130245 

52 10: DATA  '  v'  ,  290 
00, 40.582095 
, -27.520152, 
•  2",  47000,  42 
.0336515,  -3. 
1215512 

5220:DATA  "8970", 
3, 39.510754, 
87.291214,  "W 
",  1  1000,  30.  5 
93874, 85. 100 
9305 


:>230:DATA 


00,  42.-*. 


>240, 


v'  ,  280 
425060 

3,  76.4333862 
,  »Y' , 44000,  ^ 

8. 3649844,  94 
. 3318469 
DATA  "9940' , 
3. 39.3306621 
, 1 18.495637, 
'  W,  1  1000,  4  7 
.034793,  1  19. 
443853 

DATA  '  X",  270 
00,  38. 46599, 
1 22. 294  4529- 
' v : . 40000,  3- 
.  19:818,  114. 
481  "'435 
5270: DATA  '  9960'  , 
4. 42. 4250603 
.  76. 4933362, 

'  i>2  ,  :  1000,  4c 

. 4827199,  67. 
5537713 
5280: DATA  2X'  ,  250 


76 


5300; 


00, 4: . 15] 193 

. 59.533909, 

^  , 3  3002,  34. 

0 346081,  7  7. : 

446654 

PAT  A  '  7",  54  0 

00, 39. 510754 

,  87.2312:4 

DATA  ' 9970' , 

4, 24.480359? 

,  -141  .  193030 

3,  '  '^'  ,  1  1000, 

24. 1707888, - 

153.5853232 
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PC-1500  (TRS-80  PC-2)  PROGRAM  LISTING  (cont.) 


5310:DATA  "X",  300 
00, 42.443710 
4, -143.43092 
45,  "Y",  55000 
, 26.3624975, 
-128.0856445 

5320:DATA  "Z",  750 
00, 9.3245789 
,.  -138.095497 

5330:DATA  ,i9990", 
3,  57.091265, 
170. 1506789, 
"X", 1 1000, 52 
. 494404,  -173 
. 1048974 

5340:DATA  "Y", 290 
00, 65. 144030 
6, 166.531255 
,  "Z",  43000,  5 
7.262021,  152 
.2211225 
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IV.  THE  RADIO  SHACK  TRS-80  MODEL  100. 

TRS-80  MODEL  100  USER  INSTRUCTIONS. 

The  options  are  labeled  1  through  9  or  A.  Option  6  (TWO  POINT 
HEAD  &  DIST)  is  a  stand-alone  option  which  may  be  used  by  itself. 
For  LORAN  fixing,  Option  8  must  be  selected  first  to  load  the 
station  GRI's.  If  Option  8  is  not  used,  then  any  option  requiring 
the  GRI's  will  automatically  enter  Option  8  first.  Option  1  (LORAN 
FIX)  must  be  selected  at  least  once  before  Option  2  (ALTERNATE 
SOLN)  can  be  selected.  Once  the  GRI's  are  loaded,  any  of  the 
remaining  options  may  be  used  at  any  time.  When  first  run,  the  ASF 
Correction  Factors  are  set  to  zero;  if  they  are  subsequently 
entered  as  non-zero  values  they  will  remain  set  until  either 
changed  manually,  by  selection  of  Option  9,  or  until  the  program 
is  RUN  once  more.  In  all  cases,  South  latitudes  and  East 
longitudes  must  be  negative  for  input  and  will  be  labeled  by  a 
minus  sign  '-'  on  output.  For  ease  of  reference  the  options  are 
listed  by  number  rather  than  by  sequence  of  selection. 


INSTRUCTION 


DISPLAY 


INPUT 


Run  Program 


LORAN  FIX 
ITD  for  1st  GRI 
ITD  for  2nd  GRI 
Lat  of  fix 
Long  of  fix 


Select  LORAN. BA 

OPTION? 

ITD  FOR  gri  1? 
ITD  FOR  gri  2? 
LAT  =  dd^m'ss" 
LONG  =  dcWmm'ss" 
OPTION? 


[ENTER] 


1st  ITD  [ENTER] 
2nd  ITD  [ENTER] 

[ENTER] 

[ENTER] 


ALTERNATE  SOLN 
Lat  of  fix 
Long  of  fix 


LAT  =  dd°mm'ssw 
LONG  =  ddd0mm'ss" 
OPTION? 


[ENTER] 
[ENTER] 


3.  INPUT  DEST  LAT/LONG 


4.  HEAD  &  DIST  TO  DEST 
Heading  to  dest. 
Dist.  to  dest. 

OPTION? 


DEST  LAT  dd.mmss? 

Lat  of  dest. 

[ENTER] 

DEST  LONG  ddd.mmss? 

Long  of  dest. 

[ENTER] 

OPTION? 

HEADING  =  dd°mmIss" 

4 
[ENTER] 

DISTANCE  =  mmm.ff  n.mi. 

[ENTER] 
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TRS-80  MODEL  100  USER  INSTRUCTIONS  (cont.) 


INSTRUCTION 
5.  PREDICT  ITDs 

1st  ITD 
2nd  ITD 


DISPLAY 


LAT  dd.mmss? 
LONG  ddd.mmss? 
ITD  FOR  gri  1  = 

xxxxx.xx  ys 
ITD  FOR  gri  2  = 

xxxxx.xx  ys 
Hew  Lat/Lon  or 
Beturn  to  menu 

OPTION? 


INPUT 


Lat 
Long 


[ENTER] 
[ENTER] 


[ENTER] 

[ENTER] 

N  repeats  Option  5 
R  goes  to  next  line 


TWO  POINT  HEAD  &  DIST 

ORIGIN  LAT  dd.mmss? 

ORIGIN  LONG  ddd.mmss? 

DEST  LAT  dd.mmss? 

DEST   LONG    ddd.mmss? 


Heading  to  dest. 
Dist.  to  dest. 


7.  CALIBRATE 


8.  NEW  STA.  IDs 


9.  ASF  Corrections 


HEADING  =  dd0mmlss" 
DISTANCE  =  mmm.ff  n.mi 
OPTION? 


LAT  dd.mmss? 
LONG  ddd.mmss? 
ITD  FOR  gri  1? 
ITD  FOR  gri  2? 
OPTION? 


1st  GRI? 
2nd  GRI? 
OPTION? 

ASF  FOR  gri 
ASF  FOR  gri 
OPTION? 

1? 
2? 

Lat  of  origin 

[ENTER] 
Long  of  origin 

[ENTER] 
Lat  of  dest. 

[ENTER] 
Long  of  dest. 

[ENTER] 

[ENTER] 

[ENTER] 


Lat  of  fix   [ENTER] 
Long  of  fix  [ENTER] 
1st  ITD  [ENTER] 
2nd  ITD  [ENTER] 


8 
Sta.  ID  [ENTER] 
Sta.  ID  [ENTER] 


1st 
2nd 


ASF 
ASF 


[ENTER] 
[ENTER] 


A.  STOP 


A 
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TRS-80  MODEL  100  SAMPLE  PROBLEM  1 


You  are  sailing  off  of  the  coast  of  California.  Your 
destination  is  Moss  Landing  at  about  36°48'N  and  121°47'W.  Your 
present  ITD  from  station  pair  GRI  9940W  is  16019  and  your  ITD  from 
station  pair  GRI  9940Y  is  42585.  What  is  your  current  position  and 
what  is  the  heading  and  distance  from  your  current  position  to 


Moss  Landing?  We  assume 
yet  been  loaded. 


that  the  GRI ' s  9940W  and  9940Y  have  not 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


OPTION? 

DEST  LAT   dd.mmss? 

DEST  LONG  ddd.mmss? 


36.48  [ENTER] 
121.47  [ENTER] 


Input  destination 
36°48» 
121°47 ' 


OPTION? 
1st  GRI? 
2nd  GRI? 


9940W 
9940Y 


8 

[ENTER] 
[ENTER] 


Input  station  ID's, 


OPTION? 

ITD  FOR  9940W? 
ITD  FOR  9940Y? 
LAT   =  3 9° 14 '19" 
LONG  =  115°50,52" 


1 
16019 
42585 

[ENTER] 
[ENTER] 


Compute  a  fix. 


Latitude  of  fix. 
Longitude  of  fix, 


This  fix  is  in  Nevada.  There  are  always  two 
solutions  to  the  LORAN  equations.  The  navigator 
must  select  the  correct  one. 


OPTION? 

LAT   =  35° 00' 01" 

LONG  =  125° 00' 09" 


2 

[ENTER] 
[ENTER] 


Select  alternate  soln, 
Latitude  of  fix. 
Longitude  of  fix. 


This  fix  is  in  the  Pacific.  The  navigator 
knows  it  is  the  correct  fix  because  of  his 
estimate  of  where  he  should  be. 


OPTION?  4 

HEADING  =  54° 34 '11"         [ENTER] 
DISTANCE  =  190.3  8  n.mi.     [ENTER] 


Find  heading  &  distance 
to  destination. 
Heading  to  Moss  Landing. 
Distance  to  Moss  Landing 


OPTION? 


Repeat  with  new  option 
or  press  A  to  quit. 
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TRS-80  MODEL  SAMPLE  PROBLEM  2 


How  far,  and  in  what  direction,  is  Corvallis,  Oregon 
(44°34'N,  123°16*W)  from  Cupertino,  California  (37°19'N, 
122°02'W)? 


DISPLAY  CONTENTS 

OPTION? 

ORIGIN  LAT   dd.mmss? 
ORIGIN  LONG  ddd.mmss? 
DEST  LAT   dd.mmss? 
DEST  LONG  ddd.mmss? 
HEADING  =  353° 02' 59" 
DISTANCE  =  43  8.32  n.mi 


USER  RESPONSE 

COMMENTS 

6 

Two 

point  distance 

and 

heading. 

37.19  [ENTER] 

37°19' 

122.02  [ENTER] 

122° 02' 

44.34  [ENTER] 

44°34 ' 

123.16  [ENTER] 

123°16' 

[ENTER] 

[ENTER] 

OPTION? 


Repeat  with  new  option 
or  press  A  to  quit. 


TRS-80  MODEL  SAMPLE  PROBLEM  3 


Suppose  that  you  want  to  know  what  ITD's  you  would  expect  to 
receive  at  35°  North,  125°  West  and  at  36° 27 '  North,  126° 54 •  West 
from  9940W  and  9940Y.  To  determine  these  ITD's,  proceed  as  follows 
(assume  that  the  GRI ' s  have  been  loaded  as  in  Sample  Problem  1): 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


OPTION? 

LAT    dd.mmss? 

LONG  ddd.mmss? 

ITD  FOR  9940W  =  16019.35 

ITD  FOR  9940Y  =  42584.71 

Hew  Lat/Lon  or 

.Return  to  menu? 
LAT    dd.mmss? 
LONG  ddd.mmss? 
ITD  FOR  9940W  =  15572.32 
ITD  FOR  9940Y  =  43006.15 
Hew  Lat/Lon  or 

Return  to  menu? 


35  [ENTER] 
125  [ENTER] 
[ENTER] 
[ENTER] 

N 
36.27  [ENTER] 
126.54  [ENTER] 
[ENTER] 
[ENTER] 

R 


Predict 
35° 
125° 


ITD's, 


36°27 ' 
126°54' 


OPTION? 


Repeat  with  new  option 
or  press  A  to  quit. 
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TRS-80  MODEL  100  SAMPLE  PROBLEM  4 


Suppose  you  are  receiving  an  ITD  of  16308  from  9940W  and 
42800  from  9940Y.  These  ITD's  would  tell  you  that  your  location  is 
36°47,55"N  and  121°47'llnW.  However,  you  know  that  your  position 
is  benchmarked  to  be  at  36°47'36"N  and  121° 46 ' 58"W,  and  you  wish 
to  calibrate  your  Model  100  so  that  the  ITD's  of  16308  and  42800 
will  give  you  the  latter  fix  instead  of  the  former.  Proceed  as 
follows  (assume  the  GRI ' s  have  been  loaded  as  in  Sample  Problem 
1): 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


OPTION? 

LAT  dd.mmss? 
LONG  ddd.mmss? 
ITD  FOR  9940W? 
ITD  FOR  9940Y? 


36.4736 
121.4658 
16308 
42800 


[ENTER] 
[ENTER] 
[ENTER] 
[ENTER] 


Calibrate. 

36° 47  '36" 

121045I58" 


The  station  baselines  have  now  been  modified  to 
achieve  calibration.  Now,  reenter  16038  and  42800  into 
the  fixing  routine  to  test  the  calibration.  Calibration 
does  not  modify  the  permanent  station  data  stored  in  the 
Model  100.  The  calibration  remains  effective  until  the 
GRI's  are  reloaded  by  Option  8.  NOTE:  The  affect  of 
calibration  on  the  accuracy  of  fixes  far  removed  from 
the  benchmark  has  not  been  studied.  The  user  should  use 
the  calibration  option  with  caution. 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


OPTION? 

ITD  FOR  9940W? 
ITD  FOR  9940Y? 
LAT   =   36° 47' 35" 
LONG  =  121°46,59" 


16308  [ENTER] 

42800  [ENTER] 

[ENTER] 

[ENTER] 


Compute  a  fix 


Lat  of  fix. 
Long  of  fix. 


The  small  discrepancy  between  this  fix  and  the 
benchmark  is  due  to  assumptions  in  the  fixing  algorithm. 
If  higher  precision  is  required,  the  ASF  Correction 
factor  may  be  used  (see  Sample  Problem  5).  Note  that  you 
may  need  to  use  the  Alternate  Solution  Option  2  if  the 
fix  obtained  above  is  near  38052'57"N  and  116°51'11'W. 


DISPLAY  CONTENTS 
OPTION? 


USER  RESPONSE 


COMMENTS 

Repeat  with  new  option 
or  press  A  to  quit. 
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TRS-80  MODEL  100  SAMPLE  PROBLEM  5 


Suppose  you  are  receiving  an  ITD  of  12153.31  microseconds 
from  9960W  and  44451.83  microseconds  from  9960Y.  These  ITD's  would 
tell  you  that  your  location  is  44°15,05"N  and  67025'23nW.  Entering 
the  ASF  Correction  Tables  [Ref.  5,6]  with  these  readings,  the  ASF 
Correction  for  GRI  9960W  at  44°15'N  and  67°25'W  is  1.5  micro- 
seconds. Similiarly,  the  ASF  Correction  for  GRI  9960Y  at  these 

microseconds.  Enter  these  corrections  into  your 


Model  100  and  recompute  your  position. 
9960Y  using  Option  8)  . 


(Load  the  GRI ' s  9960W  and 


DISPLAY  CONTENTS 


USER  RESPONSE 


COMMENTS 


OPTION? 

ASF  FOR  9960W? 

ASF  FOR  9960Y? 

OPTION? 

ITD  FOR  9960W? 

ITD  FOR  9940Y? 

LAT   =   44°  15' 26" 

LONG  =   67° 26' 26" 

OPTION? 


9 

1.5  [ENTER] 

2.7  [ENTER] 

1 

12153.31  [ENTER] 

44451.83  [ENTER] 

[ENTER] 

[ENTER] 


ASF  Correction 


Compute  a  fix. 


Repeat  with  new  option 
or  press  A  to  quit. 


These  ASF  Corrections  have  been  entered  into  the 
computer.  The  range  of  their  validity  must  be  checked 
with  the  ASF  Correction  tables  and  reentered  if  they 
change.  NOTE:  If  you  are  using  the  ASF  Correction  Tables 
and  if  you  wish  to  calibrate  to  a  benchmark  position, 
you  should  first  enter  the  ASF  Corrections  for  your 
benchmark  and  then  use  the  calibration  Option  7  (see 
Sample  Problem  4) . 


-34- 


TRS-80  MODEL  100  PROGRAM  LISTING 

10  REM  LORAN-C  ROUTINE.  REV  07-15-83.  1530  HOURS 

20  N9=15:REM  N9  =  NO.  OF  MASTER  STATIONS 

30  C0=l/298.26 :A0=-1:PI=4.*ATN(1) :TP=PI+PI :RD=PI/180 :G$ (1) 

="" :GOTO120 

32  F0RC=1T09:C$=INKEY$:NEXTC 

33  PRINT:PRINT"Press  any  key  to  continue." 

34  C$=INKEY$:IFC$=""GOT034 

35  RETURN 

40  X=INT(M0*X+.5)/M0:RETURN 

50  P=129.04  398/T-.407  58+.64576543  8E-3*T:RETURN 

60  P=2. 7 412 979/T-. 3277 4624 -3 *T: RETURN 

70  C$=INKEY$:IFC$=""GOTO70 

80  C=ASC(C$) :C$=CHR$(C+32*(C>90)) 

90  FORC=1TOLEN(C0$) : IFC$=MID$ (C0$, C, 1 ) THENRETURN 

100  NEXTC.-GOTO70 

120  CLS: PRINT" 1-LORAN  FIX         6-TWO  PT  HEAD  &  DIST" 

130  PRINT"2-ALT.  SOLUTION     7-CALIBRATE" 

140  PRINT"3-DEST  LAT/LON      8-NEW  STA.  IDs  " 

150  PRINT"4-HEAD  &  DIST      9-ASF  Correction" 

160  PRINT"5-PREDICT  ITDs      A-STOP 

170  PRINT: PRINT"OPTION?" :C0$="1234567 89A" :GOSUB70 

180 

CLS:ONCGQSUB910, 950,1010,1040, 1090 ,  990 , 1200,190, 132 0,530 :GOTO120 

185  IFG1$<>" "THENRETURN 

190  F0=l:N=l:CLS:INPUT"lst  GRI";I$ 

200  GOSUB380:GOSUB390:IFF0=2GOTO190 

210  F0=l:N=2:INPUT"2nd  GRI " ; 1$ :GOSUB3 80 

220  IFA$=I9$THENGOSUB470 : IFF0=1GOTO250 

240  GOSUB390:IFF0=2GOTO210 

250  F1=1:F2=1:IFG1$=G2$GOTO300 

26 0  IFP ( 1 , 1 ) =P ( 1 , 2 ) ANDL ( 1 , 1 ) +L ( 1 , 2 ) GOTO310 

270  F2=-1:IFP(1,1)=P(2,2)ANDL(1,1)+L(2,2)GOTO310 

280  Fl=-1 : IFP ( 2 , 1 ) =P ( 2 , 2 ) ANDL ( 2 , 1 ) +L ( 2 , 2 ) GOTO310 

290  F2=2:IFP(2,1)=P(1,2)ANDL(2,1)+L(1,2)GOTO310 

300  PRINT"ERROR:  NO  TRIPLET. " :GOSUB32 :GOTO190 

310  PRINT:PRINT"Working":FORI=lT02:Pl=P(l,I) :Ll=L(l,I) :P2=P(2,I) 

:L2=L(2,I) 

320  GOSUB330  ;NEXTI  :G  (1)  =0:G  (  2)  =0  .-RETURN 

330  GOSUB650:B(I)=S0:Z(1,I)=Z1:Z(2,I)=Z2 

340  T=21282.3593*S0 

350  IFT>=537THENGOSUB50 

360  IFT<537THENGOSUB60 

370  T (I )=T+P: RETURN 

380    L=LEN(I$) :A$=LEFT$(I$,L-1) :B$=RIGHT$ ( 1$, 1 ) :B=ASC(B$) :B$=CHR$( 

B+32*(B>90) ) :RETURN 

3  90    N8=l: RESTORE 

400    READI9$,N0,P0,L0 

410    FORI=1TON0:READW$(I) ,D9(I) , Pi (I) , Ll (I) :NEXTI 

420    IFI9$OA$GOTO450 

430    GOSUB470:IFF0=2GOTO460 

440    RETURN 


-35- 


TRS-80  MODEL  100  PROGRAM  LISTING  (cont.) 

450  N8=N8+1:IFN8<=N9GOTO400 

460  PRINTI$;n  IS  NOT  CATALOGED" :F0=2 :GOSUB32 : RETURN 

470  FORI=1TON0:IFW$(I)=B$GOTO490 

480  NEXTI:F0=2: RETURN 

490  IFN=1THENG1$=I$:G$(1)=I$ 

500  IFN=2THENG2$=I$:G$(2)=I$ 

510  D(N)=D9(I) :X=P0:GOSUB1370:P(1,N)=X:X=L0:GOSUB1360:L(1,N)=V*RD 

520  X=P1(I) :GOSUB1370:P(2,N)=X:X=L1(I) :GOSUB1360 :L (2,N) =V* 

RD: RETURN 

530  END 

540  REM  DIRECT  SOLUTION 

550  Z8=SIN(Z1) :Z9=COS(Zl) :P8=SIN(Pl) :P9=COS(Pl) :M=-Z 8*P9 : C1=C0* 

M:C2=C0*(l-M*M)/4 

560  D=(1-C2)*(1-C2-C1*M) : P=C2* (l+Cl*M/2)/D :N=P9*Z9 

570  YY=N:XX=P8:GOSUB126  0:S1=AN:D0=S0/D:U=2*(S1-D0) :W=l-2*P*COS (U) 

:V=COS(U+D0) 

580  X=C2*C2*SIN(D0) *COS (D0* ( 2*V*V-1) ) : Y=2*P*V*W*SIN (D0) 

590  S2=D0+X-Y:S8=SIN(S2) :S9=COS(S2) : K=SQR(M*M+ (N*S9-P8*S8) ~2) 

600  P2=ATN( (P8*S9+N*S8)/K) 

610  YY=-S8*Z8:XX=P9*S9-P8*S8*Z9:GOSUB1260:S3=AN 

620  H=C1*(1-C2) *S2-C1*C2*S8*C0S(S1+S1-S2) :L2=L1+S3-H 

630  YY=-M:XX=-(N*S9-P8*S8) :GOSUB1260 : Z2=AN:RETURN 

640  REM  REVERSE  SOLN 

650  L3=L2-Ll:P3=(P2-Pl)/2:P4=(Pl+P2)/2 

660  P6=SIN(P3) :P7=COS(P3) :P8=SIN(P4) :P9=COS(P4) 

67  0    H=P7*P7-P8*P8:L=P6*P6+H*SIN(L3/2)"2:XX=SQR(L) :GOSUB1340 : D0=2* 

AN 

6  80    U=2*P8*P8*P7*P7/(1-L) :V=2*P6*P6*P9*P9/L:X=U+V: Y=U-V:T=D0/SIN ( 
D0) :D=4*T*T 

690    E=2*COS(D0) :A=D*E:C=T-(A-E)/2:N1=X*(A+C*X) :B=D+D :N2=Y* (B+E*Y) 

:N3=D*X*Y 

700    D2=C0*C0*(N1-N2+N3)/6  4:D1=C0*(T*X-Y)/4:S0=(T-D1+D2)*SIN(D0) 

710    M=32*T-(20*T-A) *X-(B+4) *Y 

720    F=Y+Y-E*(4-X) :G=C0* (T/2+C0*M/6 4 ) :Q=-F*G*TAN (L3 ) /4 

730    L4=(L3+Q)/2:L8=SIN(L4) :L9=COS(L4) 

740    YY=P6*L9:XX=P9*L8:GOSUB1260:T1=AN:YY=-P7*L9:XX=P8* 

L8:GOSUB1260:T2=AN 

750    M0=TP:X=T1+T2:GOSUB1250:Z1=X:X=T1-T2:GOSUB1250:Z2=X:RETURN 

760    REM    FIXING    ROUTINE 

770    A1=F1*SIN(A(1) ) :Bl=COS(A(l))-COS(B(l) ) :C1=SIN (B(l) ) 

7  80    A2=F2*SIN(A(2) ) :B2=COS(A(2) )-COS(B(2) ) : C2=SIN (B( 2) ) 
790    E1=Z(1,1) :IFF1=-1THENE1=Z(2,1) 

800    E2=Z(1,2) :IFF2=-1THENE2=Z(2, 2) 

810    C=B1*C2*C0S(E2)-B2*C1*C0S(E1) :S=B1*C2*SIN ( E2 ) -B2*C1*SIN (El ) 

820    K=B2*A1-B1*A2:R=SQR(C*C+S*S) : YY=S:XX=C:GOSUB1260 :G=AN 

830    XX=K/R:GOSUB1350:Z=G+A0*AN:YY=B2:XX=C2*COS(Z-E2)+ 

A2:GOSUB1260:S0=AN 

840    IFF2=1THENP1=P(1, 2) :Ll=L(l,2) 

850    IFF2=-1THENP1=P(2,2) :Ll=L(2r2) 

860    Z1=Z:GOSUB550:P0=P2:L0=L2:P=ATN(TAN(P0)/(1-C0)) 
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TRS-80  MODEL  100  PROGRAM  LISTING  (cont.) 

870  P=P/RD:X=L0/RD:M0=360:GOSUBl250:L=X:IFL>180THENL=L-360 

880  X=P:GOSUB1280:P$=C$:X=L:GOSUB1280 

890  PRINT: PRINT "LAT   =  " ;P$ :PRINT"LONG  =  " ;C$ : RETURN 

900  REM  OPT  1 

910  GOSUBl 85 :CLS: FORI =1T02 

920  PRINT" ITD  FOR  " ;G$ (I) ; : INPUTA:A=A+G ( I) -D ( I ) -T (I) : IFABS (A)<T ( I ) 

GOTO940 

930  PRINT"ITD  NOT  VALID  FOR  " ;G$  (I)  :GOSUB32 :GOTO920 

940  A(I) =A/21295 . 8736 :NEXTI :GOSUB770 :GOSUB32 :RETURN 

945  REM  OPT  2 

950  IFA0=1THENA0=-1:GOTO970 

960  A0=1 

970    GOSUB770:GOSUB32:RETURN 

980    REM    OPT    6 

990  CLS:INPUTnORIGIN  LAT    dd.mmss "; P: INPUT"ORIGIN  LONG 

ddd.mmss" ;L 

1000  X=P:GOSUB1370:P0=X:X=L:GOSUB1360 :L0=V* 

RD:GOSUB1010 : GOSUBl 04 0 : RETURN 

1010  CLS:INPUT"DEST  LAT    dd.mmss" ; P: INPUT"DEST  LONG  ddd.mmssn;L 

1020  X=P:GOSUB1370 :P5=X:X=L:GOSUB1360 :L5=V*RD:RETURN 

1030  REM  OPT  4 

104  0    Pl=P0 :L1=L0 : P2=P5 :L2=L5 :GOSUB650 

1050    M0=360:X=Z1/RD:GOSUB1250:GOSUB1280 

1060    CLS:PRINT"HEADING    =    " ; C$ :M0=100 : X=3443 . 917  * 

S0:GOSUB40:PRINT"DISTANCE   =    " ;X;n    n.    mi." 

1070    GOSUB32: RETURN 

10  80    REM   OPT    5 

1090    GOSUBl 85 :GOSUB117  0: PRINT :F0RK=1T02:I9=T(K)+D(K) 

1100    P2=P(2,K) :L2=L(2,K) :GOSUB650 : 1=3 :GOSUB340 : I9=I9+T (3) 

1110    P2=P(1,K) :L2=L(1,K) :GOSUB650 : 1=3 :GOSUB340 : I9=I9-T(3) 

1120    M0=100:X=I9:GOSUB40 

1130  PRINT"ITD  FOR  ";G$(K);n  =  ";X:NEXTK 

1140  PRINT: PRINT"New  lat/long  or  Return  to  menu?" :C0$="NR" :GOSUB70 

1150  IFC=1GOTO1090 

1160  RETURN 

117  0    CLS:INPUT"LAT        dd.mmss" ; PI : INPUT"LONG   ddd.mmss" ;Ll 

11 80    X=P1 :GOSUB1370 : P1=X :X=L1 :GOSUB1360 :Ll=V*RD:RETURN 

1190    REM   OPT   7 

1200  GOSUBl 85 :GOSUB1170 : FORK=lT02 : I$=G1 $ : IFK=2THENI$=G2 $ 

1210  CLS:PRINT"ITD  FOR  ";I$;:INPUT  19 : I9=I9-D(K) 

1220  P2=P(2,K) :L2=L(2,K) :GOSUB650 : 1=3 :GOSUB340 : I9=I9-T (3) 

1230  P2=P(1,K) :L2=L(1,K) :GOSUB650 : 1=3 :GOSUB340 :T (K) =I9+T (3) +G (K) 

1240  NEXTK:RETURN 

1250  X=X-M0*INT(X/M0) :REM  MOD  FCTN 

1260  AN=ATN(YY/(XX~1E-9*(XX=0) ) )-PI*(XX<0) :RETURN:REM  QATN 

1270  X=ATN( (1-C0)*TAN(X) ): RETURN 

1280  C$="  ":IFX<0THENC$=n-":X=-X 

1290  X=X+1/7200:X0=INT(X) :C$=C$+STR$(X0)+"  " 

1300  X=60*(X-X0) :X0=INT(X) :X$=STR$ (100+X0) : C$=C$+RIGHT$ (X$f 2) +"  " 
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TRS-80  MODEL  100  PROGRAM  LISTING  (cont.) 

1310  X=60*(X-X0) :X0=INT(X) :X$=STR$ (100+X0) : C$=C$+RIGHT$ (X$, 2) +  " 

": RETURN 

1320  GOSUB1 85 : CLS : PRINT "ASF  FOR  " ;G$ (1) ;: INPUTG (1) 

1330  PRINT"ASF  FOR  " ;G$ (2) ; : INPUTG ( 2) rRETURN 

1340  ZZ=SQR(1-XX*XX) :AN=ATN (XX/ ( ZZ-1E-9* ( ZZ=0) ) ) : RETURN: REM  ASIN 

1350  AN=ATN(SQR(1-XX*XX)/(XX-1E-9*(XX=0) ) ) -PI*(XX<0) :RETURN:REM 

ACOS 

1360  S=SGN(X) :X=ABS(X) :  H=INT (X) :M0=1 :GOSUB1250 :V=X* 

100:X=V:GOSUB1250 

1365  V=S*( (100*X/60+INT(V) J/60+H) :RETURN 

1370  GOSUB1360 :X=V*RD:GOSUB127  0 : RETURN 

5000  DATA 

"4990 ", 2, 16. 4443 95 , 16 9. 30312 ,nX", 11000 ,20. 144916 ,155. 53 097 

5010  DATA  "Y", 29000, 28. 234177, 178. 17302 

5020  DATA 

"5930 ' ', 2, 46. 482 71 99, 67. 553 771 3,"X", 11000, 41. 1511 93, 6 9. 583 90 9 

5030  DATA  "Y", 25000, 46. 46463218, 53. 102816 

5040  DATA  "5970 ",3,36. 11057 97, -129. 2027 27  9, "W" , 11000, 42. 4437104,- 

143 .  4309245 

5050  DATA  "X ",31000, 35. 0223 871, -126. 32267 4,MZ", 42000, 263624975, - 

128.0856445 

506  0  DATA 

"5990  ",3,51. 57587 8, 122. 220224, "X", 11000, 55. 26 20 851, 131. 151 9648 

507  0  DATA 

"Y", 27 000, 47. 0347 99, 11 9. 443 953, "Z", 41 000, 50. 3629731, 127. 2129043 

50  80  DATA 

"7 930 ",3,59. 5917 27, 45. 1027 47, nW", 11000, 6 4. 5426 58, 23. 55217 5 

5090  DATA 

"X ",21000, 62. 17 5964, 7. 0426 53 8, "Z", 43 000, 46. 46321 8, 53. 102 816 

5100  DATA  " *7 930 ",3,24. 17 07 888, -153. 5853232, "X", 11000, 42.  4437104, - 

143.4309245 

5110  DATA  " Y ", 30 000, 26. 3624 975, -12 8. 0  856 445," Z", 490 00, 9. 32457 89,- 

138.095497 

5120  DATA 

"7 960 ",2,63. 1942 814, 142. 48319, "X", 11000, 57. 26 2021, 152. 2211225 

5130  DATA  "Y", 26000, 55. 2620851, 131. 1519648 

5140  DATA  "7970", 4, 62. 175964, 7. 0426538, "W", 26000, 54. 4829872, - 

8.1736312 

5150  DATA  "X", 11000, 68. 380615, - 

14.27  47, "Y", 460 00, 6 4. 5426 5 8, 23. 55217 5 

5160  DATA  "Z", 60000, 70. 545261, 8.435869 

517  0  DATA 

"7 980 ",4,30. 593 874, 85. 1009305, "W", 11000, 30. 433301 8, 90. 49436 

5180  DATA 

"X ",23 000, 26. 3155006, 97. 50000 93, "Y", 43 0  00, 27. 015 83 93, 80. 06 53 429 

5190  DATA  "Z", 59000, 34. 0346081, 77. 5446654 

5200  DATA  "7 990 ", 3, 3 8. 52205 87, -16. 4306159, "X", 11000, 35. 31207 87,- 

12.3130245 
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TRS-80  MODEL  100  PROGRAM  LISTING  (cont.) 

5210  DATA  "Y ", 29000 , 40. 582095 , -27. 520152 ,  "Z", 47000, 42. 0336515 ,- 

3.1215512 

5220  DATA 

"8970", 3, 39.5107 54,87. 291214, "W", 11000 , 30 . 593  874 , 85 .100  9305 

5230  DATA 

"X ",28000, 42. 42506 03, 76. 4933 862, "Y", 44000, 48. 36 49844, 94. 331 846 9 

5240  DATA 

"9940 ",3,3 9. 3306621, 11 8. 495637, "W", 11000, 47. 0347 99, 11 9. 443 9b3 

5250  DATA 

"X",27  000,3  8.46  56  99,122.2944529,"Y",40000,35.191818,114.4817  435 

527  0  DATA 

"9960 ",4,42. 42506 03, 76. 4933 862, "W", 11000, 46. 48271 99, 67. 5537713 

52  80  DATA 

"X ",25000, 41. 1511 93, 6 9. 583909, "Y", 3 9000, 34. 03460 81, 77. 5446654 

5290  DATA  "Z " , 54000 , 39 . 510754 , 87 . 291214 

5300  DATA  " 9970 ",4,24. 4803597, -141. 1930303, "W", 11000, 24. 17 07 888,- 

153.5853232 

5310  DATA  "X", 30000, 42. 44371 04, -143. 43 09245, "Y", 55000, 26. 3624975, ■ 

128.0856445 

5320  DATA  "Z " ,75000 , 9 . 32457 89 ,-13 8. 095497 

5330  DATA  "9990 " ,3 , 57 . 09126  5 ,17  0 .15067 89 , "X" , 11000 , 52  .  494404  ,- 

173.1048974 

5340  DATA 

"Y", 29000, 6 5. 1440306, 166. 531255, "Z", 43000, 57. 26 2021, 152. 2211225 
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V.  MEMORY  REQUIREMENTS 


HP-75C 

PC 
(TRS- 

:-1500 
•80    PC- 

•2) 

TRS- 80 
MODEL    100 

Memory    (bytes) 

for : 
Program 

7747 

7189 

7220 

Variables 

4920 

? 

47  86 

TOTAL 

12667 

9 

• 

12006 
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DMA 
Stock  Number  Title 

LCPUB2211100C  East  Coast  Canada  (Pairs  5930-X,Y) 

LCPUB2211200C  Northeast  USA  (Pairs  9960-W,X, Y, Z) 

LCPUB2211300C  Great  Lakes  (Pairs  8970-X,Y,Z) 

LCPUB2211400C  Southeast  USA  (Pairs  7980-W,X, Y, Z) 

LCPUB2212100C  West  Coast  of  USA  (Pairs  9940-W,X,Y) 

LCPUB2212200C  West  Coast  Canada  (Pairs  5990-X,Y,Z) 

LCPUB2212300C  Gulf  of  Alaska  (Pairs  7960-X,Y) 

LCPUB2212400C  North  Pacific  (Pairs  9990-X,Y,Z) 
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Appendix:   LORAN-C  STATION  DATA 

The  following  list  contains  the  pertinent  parameters  for  each 
Loran-C  station  pair.  This  list  was  compiled  from  the  data  in 
Reference  4.  Each  column  contains  the  following  information: 

1.  The  Loran-C  station  pair  designator. 

2.  The  coding  delay. 

3.  The  master  station  latitude. 

4.  The  master  station  longitude. 

5.  The  slave  station  latitude. 

6.  The  slave  station  longitude. 

In  this  list,  positive  longitudes  are  West,  negative 
longitudes  are  East,  positive  latitudes  are  North  and  negative 
latitudes  are  South.  In  columns  3  through  6  the  latitudes  and 
longitudes  appear  to  be  in  decimal  form,  but  the  actual  format  is 
DDD.MMSSFFF  where 

DDD  designates  degrees, 
MM  designates  minutes, 
SS  designates  seconds,  and 
FFF  designates  thousandths  of  seconds. 
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LORAN-C    STATIONS 


GEL  ££  MS    LAT  MS    LQN  SS    LAT  SS    LQN 

4990X  11000  16.4443950         169.3031200  20.1449160  155.5309700 

4990Y  29000  16.4443950         169.3031200  28.2341770  178.1730200 

5930X  11000  46.4827199         067.5537713  41.1511930  069.5839090 

5930Y  25000  46.4827199         067.5537713  46.4632180  053.1028160 

5970W  11000  36.1105797  -129.2027279  42.4437104  -143.4309245 

5970X  31000  36.1105797  -129.2027279  35.0223871  -126.3226741 

5970Z  42000  36.1105797  -129.2027279  26.3624975  -128.0856445 

5990X  11000  51.5758780         122.2202240  55.2620851  131.1519648 

5990Y  27000  51.5758780         122.2202240  47.0347990  119.4439530 

5990Z  41000  51.5758780         122.2202240  50.3629731  127.2129043 

7930W  11000  59.5917270         045.1027470  64.5426580  023.5521750 

7930X  21000  59.5917270         045.1027470  62.1759640  007.0426538 

7930Z  43000  59.5917270         045.1027470  46.4632180  053.1028160 

*7930X  11000  24.1707888  -153.5853232  42.4437104  -143.4309245 

*7930Y  30000  24.1707888  -153.5853232  26.3624975  -128.0856445 

*7930Z  49000  24.1707888  -153.5853232  09.3245789  -138.0954970 

7960X  11000  63.1942814         142.4831900  57.2620210  152.2211225 

7960Y  26000  63.1942814         142.4831900  55.2620851  131.1519648 

7970W  26000  62.1759640  +007.0426538  54.4829872  -008.1736312 

7970X  11000  62.1759640  +007.0426538  68.3806150  -014.2747000 

7970Y  46000  62.1759640  +007.0426538  64.5426580  +023.5521750 

7970Z  60000  62.1759640  +007.0426538  70.5452610  +008.4358690 

7980W  11000  30.5938740         085.1009305  30.4333018  090.4943600 

7980X  23000  30.5938740         085.1009305  26.3155006  097.5000093 

7980Y  43000  30.5938740         085.1009305  27.0158393  080.0653429 

7980Z  59000  30.5938740         085.1009305  34.0346081  077.5446654 

7990X  11000  38.5220587  -016.4306159  35.3120787  -012.3130245 

7990Y  29000  38.5220587  -016.4306159  40.5820950  -027.5201520 

7990Z  47000  38.5220587  -016.4306159  42.0336515  -003.1215512 

8970W  11000  39.5107540         087.2912140  30.5938740  085.1009305 

8970X  28000  39.5107540         087.2912140  42.4250603  076.4933862 

8970Y  44000  39.5107540         087.2912140  48.3649844  094.3318469 

9940W  11000  39.3306621         118,4956370  47.0347990  119.4439530 

9940X  27000  39.3306621         118.4956370  38.4656990  122.2944529 

9940Y  40000  39.3306621         118.4956370  35.1918180  114.4817435 

9960W  11000  42.4250603         076.4933862  46.4827199  067.5537713 

9960X  25000  42.4250603         076.4933862  41.1511930  069.5839090 

9960Y  39000  42.4250603         076.4933862  34.0346081  077.5446654 

9960Z  54000  42.4250603         076.4933862  39.5107540  087.2912140 

9970W  11000  24.4803597  -141.1930303  24.1707888  -153.5853232 

9970X  30000  24.4803597  -141.1930303  42.4437104  -143.4309245 

9970Y  55000  24.4803597  -141.1930303  26.3624975  -128.0856445 

9970Z  75000  24.4803597  -141.1930303  09.3245789  -138.0954970 

9990X  11000  57.0912265  +170.1506789  52.4944040  -173.1048974 

9990Y  29000  57.0912265  +170.1506789  65.1440306  +166.5312550 

9990Z  43000  57.0912265  +170.1506789  57.2620210  +152.2211225 
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